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The AU Microscopii Debris Disk: Multiwavelength Imaging and 
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James R. Graham^ 

ABSTRACT 

Debris disks around main sequence stars are produced by the erosion and 
evaporation of unseen parent bodies, which are potential building blocks of plan- 
ets. Such planets may imprint dynamical signatures on the structure of the debris 
disk. AU Microscopii (GJ 803) is a compelling object to study in the context of 
disk evolution across different spectral types, as it is an M dwarf whose near 
edge-on disk may be directly compared to that of its A5V sibling (3 Pic. We re- 
solve the disk from 8-60 AU in the near-IR JHK' bands at high resolution with 
the Keck II telescope and adaptive optics, and develop a novel data reduction 
technique for the removal of the stellar point spread function. The point source 
detection sensitivity in the disk midplane is more than a magnitude less sensitive 
than regions away from the disk for some radii. We measure a blue color across 
the near-IR bands, and confirm the presence of substructure in the inner disk. 
Some of the structural features exhibit wavelength-dependent positions. 

The disk architecture and characteristics of grain composition are inferred 
through modeling. Previous efforts have modeled the dust distribution through 
a variety of means, ranging from power law models to calculations of steady- 
state grain dynamics. Recent measurements of the polarization properties of the 
scattered light indicate the presence of porous grains. The scattering proper- 
ties of these porous grains have a strong effect on the inferred structure of the 
disk relative to the majority of previously modeled grain types. We approach 
the modeling of the dust distribution in a manner that complements previous 
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work. Using a Monte Carlo radiative transfer code, we compare a relatively 
simple model of the distribution of porous grains to a broad data set, simulta- 
neously fitting to midplane surface brightness profiles and the spectral energy 
distribution. Our model confirms that the large-scale architecture of the disk is 
consistent with detailed models of steady-state grain dynamics. Here, a belt of 
parent bodies from 35-40 AU is responsible for producing dust that is then swept 
outward by the stellar wind and radiation pressures. We infer the presence of 
very small grains in the outer region, down to sizes of ~0.05 /xm. These sizes are 
consistent with stellar mass- loss rates M^, <^ 10^ M©. 

Subject headings: circumstellar matter — planetary systems: protoplanetary 
disks — stars: individual (AU Mic) — stars: low-mass, brown dwarfs — instru- 
mentation: adaptive optics — techniques: image processing 



INTRODUCTION 



In the past two decades, observations enabled by new technology have substantially 
broadened our understanding of dust surrounding main-sequence stars. These tenuous 
disks of dust, found aro und nearby stars, are fuiidamentally linked to the processes of 



star and planet for mation (IBackman fc Parescdll993l : lLagrange et al.ll200Cll : IZuckermarul2001 



Meyer et al.l 120071 ). Current scenarios favor casting "debris disks" in the later stages of the 
systems' formati on, after most of the pr i mordial gas and du st have dissipated (on timescales 
of 1-10 Myr; e.g. IZuckerman et al.lll995l : iHaisch et al.ll200ll ). The debris is freshly nourished 
by the sublimatio n, evaporation, an d collisional destruction of orbiting parent bodies (for a 
recent review, cf. iMeyer et al.l 120071 ). Since a portion of these parent bodies may take part 
in the accretion of cores into planets — which may in turn gravitationally perturb the dust 
— circumstellar debris disks offer evidence for the presence of planets and insight into their 
formation. 

Debris disks are usually discovered by sensing the thermal re-radiation of abs orbed 
starlight in the far IR (IBackman fc Parescelll993l : lDecin et al.ll2003l : iBryden et al.ll2006l ). The 
spectral energy distributions (SEDs) of these systems, particularly the range from the IR to 
the sub-mm, are comr nonly used to inf er the (often cold) dust temperature and extent of 
disk inner clearing (e.g. lChen et al.ll2005l ). However, in the absence of additional information 



significant degeneracies exist in modeling the dust disks with SEDs (IMoro-Martm et al.ll2005 



and references therein) . Severe ambiguities between disk morphology and the grains' physical 
characteristics are removed with spatially resolved images of the thermal and scattered light; 
however, the faintness of the dust and often-large dynamic range between the star and disk 
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present significant obstacles to their acquisition. Resolved images are both rare and valuable. 



Images at multiple wavelengths are yet more powerful. Spatially resolving the debris 
across several wavebands provid es definitive evidence for the size, space, and compositiona l 
distributions of the grains (e.g. iPantin et al.l 119971 : IWyattI l2006l : iGolimowski et al.l l2006l ). 
Coupling such data with detailed models of grain dynamics can expose the fundamental 
physical phenomena governing the production and diffusion of dust throughout the disk. 
Physical models trace the characteristics of grains (i.e. composition and geometry) to the 
underlying parent body distribution, which in turn enables inferences regarding the co- 
agulation history of solid material. These parent bodies may be corralled into a specific 
architecture through the gravity of n earby planets. Moreover, planets can imprint their 



2OOOI : iKuchner fc Holmanl l2003l : iKenyon fc Bromley 



signatures on the dust disk itse l f fe.g. iLiou fc Zooklll999l: IWvatt et al. Il999l : lOzernoy et al. 



2004 



WyattI I2OO6I ) . Tracing the mech- 



anisms of grain dynamics through resolved imaging also addresses crucial questions about 
the nature of dust production in these systems — is the dust distribution steady-state, 
or is the evolution of debris governed by transient events, such as rare collisions between 
massive bodies or extrasolar analogs to the Late Heavy Bombardment? Again, the images 
that are critical to addressing these issues are rare, though the number o f resolved systems 
has increased rapidly in the past few years (for recent tabulations, see iKalas et al.l 12006 



and lMeyer et al. 



20071 ). With ever- increasing sample sizes, we have entered an era where we 



can compare the structure and evolution of circumstellar dust with quantities which shape 
the paths of planet formation, such as stellar mass and metallicity. 

Contemporary work has revealed an exciting laboratory for the detailed study of cir- 
cumstellar debris. AU Microscopii (GJ803) is a nearby star that harbors an optically thin 



debris disk. The spectral type in the literature ranges f rom dM0-2.5Ve (jJoy fc Abtlll974 



Linsky et al.l Il982l : iKeenanI Il983l : iGliese fc Jahreisd Il995l ) . It is a member of the /3Pictoris 



moving group, and as such is one of the youngest (12]'^4Myr; iBarrado v Navascues et a 



I999I : IZuckerman et al.ll200ll ) and among the closest (9.94 ± 0.13 pc; iPerryman et al.l 119971 ) 



of the known resolved debris disks. It is especially attractive to study in the larger context 
of disk evolution as a function of stellar mass. In terms of mass, M dwarfs like AU Mic con- 
stitute the dominant stellar component of the Galaxy. Despite the abundance of such stars, 
AU Mic is currently unique among the resolved debris systems. It is a touchstone for study- 
ing the evolution of circumste llar disks around low- mass stars. AU Mic has a well-chronicled 
history of flare activity (e.g. iRobinson et al.ll200ll ). It is likely that strong stellar activity 
has a significant effect on dust dynamics and lifetimes (jPlavchan et al.ll2005l ). Notably, we 
can make direct comparisons across spectral types given AU Mic and its sibling (3 Pic, the 
archetypal A star with an edge-on debris disk. 
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The scattered light of the dust around AU Mic was discovered by iKalas et aLl (120041 ). 
who used seeing-hmited i?-band coronagraphic imaging to resolve it into an extended, near- 
edge-on disk. Additior ial h i gh-re solution st udies follo wed, with adap t ive op tics (AO) obser- 



vations at H-hand by iLiul (12004 hereafter |L0J) and iMetchev et all (12005 



hereafter 1MQ5|) 



as we ll as in the visible with the Hubble Space Telescope (HST) by lKrist et al. J2OO5I . here- 



after |K05|). These detailed images reveal a very thin midplane (FWHM~2AU) with an 
inner disk closely aligned with the line of sight {6 ~ 0.5° for r < 50 AU: lK05l ). The midplane 
surface brightness profiles show a break around 35-45 AU, with brightness decreasing more 
sharply at larger projected distances. There are slight asymmetries in the overall brightness 
of the profiles between the two disk ansae. Further, the midplane exhibits substructure 
on smaller scales, including localized enhancements and deficits in brightness, and vertical 
deviations of midplane positions away from that of a uniform disk. Comparison of these 
structures between datasets is required to confidently reject image processing artifacts. The 
origin of the small-scal e str ucture remains unexplained. A striking feature of the scattered 
light disk is its color. lK05l found a blue color in the visible with HST, along with an ap- 
parent color gradient — the disk is increasi ngly b lue from 30-60 AU. The disk is also blue 
from the visible to the H band, as noted by lM05l . This is unlike many of the disks resolved 



to da te, which are neutral or red scatterers like (3 Pic ( Golimowski et al.ll2006l: Mever et al. 



20071 ). We note, however, th at two other recently imag ed disks, HD 32297 JKalasll2005h ami 
HD 15115 (IKalas et al.l 120071 ). also scatter blue between the optical and near-IR. Such data 
are important because images at each new wavelength serve to further constrain the sizes, 
composition, and structure of grains through analysis of their scattering properties. 

While the scattered light imaging of the disk has attracted attention in the past few 
years, there has also been significant study of the grains' thermal emission in the SED. The 
cold dust around AU Mic was first identified by th e Infrared Astronomical Satel l ite (IRAS), 



as ex cess emission at 60 /im was weakly detected (IMathioudakis fc Doyld Il99ll : ISong et al. 



2OO2I ). Recently, sensi tive measureme nts of the 850 /zm flux confirmed the presence of cold 
dust around the star (ILiu et al.ll2004l ). The Syit zer Space Telesco pe has further constrained 
the thermal emission of the dust in the mid-IR (jChen et al.l 120051 ). 



Results of these studies have been taken as indirect evidence for the presence of planets 
around the star. The gravitational influence of such planets provides a possible mechanism for 
generating the substructure in the disk. In particular, mean-motion reso nances from a planet 



can trap dust to produce radial, azimuthal, and vertical stru cture (e.g. lOzernoy et al. 



2000 



Wahhaj et al.ll2003l : iThommes fc Lissaueiil2003l : IWyattll2006l ). though to date, no models of 
these mechanisms have been applied to the substructure of AU Mic. Further, the shallow 
surface brightness profiles in the inner disk and lack of thermal excess in the 10-20 /xm 
region of the SED have suggested an inner clearing of dust, perhaps maintained again by 
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the presence of an inner planet (e.g. iRoques et al.l Il994l : iMoro-Martm fc Malhotral l2005l ). 
Nonetheless, it is imperative to understand the detailed physics affecting the distribution of 
dust grains around this disk prior to establishing the presence of unseen planets. The gravity 
of such perturbers is just one of many effects which may shape the dust distribution, including 
grain-grain collisions, forces from both stellar radiation and wind, and gas drag. Because the 
disk is along the line-of-sight to the star, AU Mic is a favorable tar g et for using absorption 
spectroscopy to search for remnant gas in the disk. iRoberge et al.l (120051 ) placed limits on 
the column of H2 toward the star using the Far Ultraviolet Spectroscopic Explorer (FUSE), 
concluding that giant planets were unlikely to have formed given the rapid dissipation of 
primordial gas. Different groups have placed pla net detection limits aroun d the star in very 
narrow (r < 2" ^ 20 AU; iMasciadri et al.l |2005| ) and slightly wider fields (INeuhauser et al. 



2OO3I : IMOSI ). and to date none have been directly detected. Certainly the question of planets 
in this system warrants further study, as these different lines of evidence have not been 
resolved. 

Comprehending the physics that sculpt the dust distribution in this debris disk remains 
a key step for not only determining the existence of planets, but also for inforrn i ng us 
abo ut the evolution of solid material around stars in general. lAugereau fc BeustI (120061 ) 
and IStrubbe &: Chiang (120061 . hereafter ISCOq ) have both investigated the observed break in 
the midplane surface brightness profile and color gradient in the outer disk, and hypothesize 
that a ring of parent bodies (near the break in surface brightness profile) acts as a s ource of 
dust grains, which are subsequently swept outward by a strong stellar wind. ISC06l provide 
a detailed physical model for AU Mic for this scenario. The steady-state spatial and size 
distributions of dust grains are determined by a small set of physical processes (e.g. collisions, 
radiation, wind, etc.), and with other assumptions about grain properties, the model can 
reproduce observations of scattered light profiles and the SED. 

As illustrated by the work done to date, we have constraints on the distribution of 
dust using scattered light images coupled with the SED. However, in this edge-on system, 
a degeneracy remains between the scattering properties of grains and their supposed spatial 
distribution. The advantage of measurements in polarized light, rather than total intensity, 
is that they give complementary information in the grain scattering properties (for grain sizes 
a such that x ~ 27ra/A ~ 1) and thus reduce this degeneracy. The polarization pr o pertie s 
of AU Mic's disk have recently been measured with HST/ ACS by I Graham et al.l (120071 ). 
The disk exhibits strong gradients in linear polarization, and their modeling of the flux 
of polarized visible light indicates the presence of small porous grains and an architecture 
consistent with a significant (>300:1 in vertical optical depth to scattering) inner clearing. 



In this paper, we report and analyze multiband AO observations of the scattered light 
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disk of AU Mic in the near IR. Our high-contrast images were processed with a novel data 
reduction technique which aims to mitigate the effects of point-spread function variabihty 
on ground-based observations. We report on the observed colors and color gradients of 
the disk brig htnes s profiles for our images, and also for reprocessed HST data previously 
pr esented bv IK05. Gu ided by the recent characterization of grain scattering properties 



by lGraham et al.l (120071 ). we simultaneously fit a dust model to near IR and visible scattered 
light data as well a s the SED. With this model of dust distribution, we then check the 
consistency of lSCOa s physical model for the disk architecture with our empirical results. We 
examine the question of whether the inner disk is populated, as expected if corpuscular drag 
forces can draw grains inward before they are pulverized through collisions. We document 
our observational and data reduction techniques in ^ we present our observational results 
in 113] and analyze dust models in ^ Discussion and conclusions follow in §3 



2. OBSERVATIONS & REDUCTION 

The achievable spatial resolution is a prime motivator for imaging debris disks in the 
optical and near infrared — resolved images here highlight structural details and dust prop- 
erties. However, at these wavelengths the star overwhelms the light scattered by the cir- 
cumstellar debris. This contrast ratio characterizes a fundamental observational challenge in 
such studies of debris systems. To meet this challenge, we use AO to concentrate the star's 
light and employ a coronagraph to occult the resulting stellar image, thereby increasing 
sensitivity to the disk. Here, we detail our use of AO coronagraphy and image processing 
techniques to study the circumstellar dust ( §2.1f[2l2!) . We describe calibration of the AO 
data in §2.3[ and our use of HST images in §2.41 



2.1. Near Infrared Imaging 

We observed AU Mic on the nights of 2004 Aug 29-30 with the Keck II AO system and 
a coronagraphic imaging mode of the NIRC2 camera. AU Mic is sufficiently bright to serve 
as its own reference for adaptive wavefront correction. After wavefront compensation, the 
on-axis starlight is blocked by a focal-plane mask. Diffraction effects are then suppressed by 
a pupil-plane Lyot stop, and the light is reimaged onto a 1024 x 1024 pixel Aladdin detector. 
The narrow- field mode of the camera was used, at a scale of ^10 mas pixel"^. Exposure 
times were 30-60 s in JHK' (60 x 1 s in L'), and the filter was cycled after a few short 
exposures in each band. Exposures where the disk was aligned with a diffraction spike were 
discarded. The filters, total integration times, and focal masks are listed in Table [H We used 
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a 0'/5 radius focal plane mask on the second night, however the data in the region around 
the edge of this mask were discarded due to large residual subtraction errors. Calibration 
employed standard bias subtraction and flat fielding techniques. 

Imaging faint circumstellar material requires suppression of the stellar light, whose dis- 
tribution is given by the on- axis point spread function (PSF). While the coronagraph sup- 
presses much of the starlight, a fraction leaks through the system and must be subtracted 
from the data in post-processing. The relative success of this operation is linked to PSF 
stability. In observations through an altitude-azimuth telescope like Keck, the image plane 
rotates relative to the telescope with parallactic angle. However, features in the PSF that 
are produced by the telescope (e.g. diffraction spikes) maintain a fixed orientation relative 
to it. We disabled the instrument rotator of NIRC2 to additionally fix the orientation of any 
features arising from camera aberrations. With the disabled rotator, the PSF orientation 
is fixed relative to the detector while the edge-on disk appears to rotate around the stellar 
image with time. This allows us to disentangle the image of the disk from features in the 
diffraction pattern of the star, as well as use the stable orientation of the stellar PSF on 
the detector for more accurate subtraction. "Roll subtraction" is a general observational 
methodolo gy that has been applied with success to AO imaging of circumstellar disks (e.g. 



L04l : IMOSI ). An account of our version of the roll subtraction technique follows. 



2.2. PSF Subtraction Technique 



The time variability of the AO PSF is often the limiting factor in detecting faint cir- 
cumstellar material from the ground. Rather than subtracting the stellar PSF from non- 
contemporaneous observations of reference stars, we wi sh to directly exploi t info rmation from 



the sc ience target exposures. Techniques developed by lVeran et al.l (119971 ) and lSheehy et al. 



( 120061 ) seek to model the AO system to estimate the PSF, using telemetry data and crowded 
field imaging, respectively. However, these techniques do not yet have demonstrated appli- 
cability to high-contrast imaging. For roll subtraction, mitigation of AO PSF variability 
can poten tially improve contr ast. We note in particular the technique's advancement in this 
regard by iMarois et al.l (120061 ) for the case of point source detection. 



We developed an algorithm suited to the self-subtraction of time-varying AO PSFs and 
applied it to the reduction of our near-IR images. In essence, because the PSF remains in 
a fixed orientation while the edge-on disk rotates, one may use PSF information spanning 
several exposures to estimate the stellar PSF (effectively removing the disk). Our refinement 
to the procedure is to estimate the PSF for each exposure, rather than for the ensemble. The 
PSF estimates are then used to subtract the stellar image from each frame. The residuals are 
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transformed to a common sky orientation and then combined to estimate the star-subtracted 
object field (Fig. [1]), corresponding to the disk emission. We detail our procedure, including 
the refinements for tracking changes in the time-variable AO PSF, in Appendix \M 



2.3. Calibration 



Immediately after bias subtraction and fiat fielding, we corrected the camera's geometric 
distortion in each image. However, no attempt was made to calibrate for the scale and 
orientation of the detector on the sky, as any differences from nominal values are expected 
to be minor. When visually comparing with HST imaging ( §2.4p . we found a 0?1 rotation 
between the two sets. 

Our calibrations for photometry require measurement of zero point and PSF. We brack- 
eted our observations of AU Mic with photor netric standard s t ars SJ9170 and HP 205772 
on the first night, and GJ811.1 on the second (lElias et al.lll982l : iPersson et al.lll998l ). These 
observations were used to determine the photometric zero point in each band. The stars 
were positioned outside the coronagraphic spot, and aperture photometry was used to mea- 
sure the stellar brightness. Since large aperture radii could be used on the well-exposed 
stellar images, no encircled energy corrections were used in determining zero points. For all 
measurement s, we applied an airma ss correction assuming extinction values appropriate for 
Mauna Kea (IKrisciunas et al.l 119871 ). 



To measure the scattering properties of the dust, we express the disk fiux relative to that 
of the star. However, measuring the brightness of AU Mic is a challenge with this instrumen- 
tal configuration because it is too bright for direct unocculted imaging, while observations 
using the focal plane mask complicate calibration. Here we adopt 2MASS ph otometry for the 
stella r brightness (J = 5.436 ± 0.017, H = 4.831 ± 0.016, K = 4.529 ±0.020; Iskrutskie et aL 
20061 ). We ignore the small (~ few centi-mag) color correction when transforming from 



K to K', as we expect disk photometry to be dominated by other errors. To determine 
the brightness of the star in L', we first measured the brightness of GJ 811.1 through the 
partially-transmissive focal plane spot in L. We then applied a color correction to compute 
the zero point (behind the spot) in L' , and measured the brightness of the AU Mic using 
the same photometric aperture, finding 4.38 ±0.04 mag in L'. We note that our attempts 
to express the relative disk and stellar fiux ratio should take AU Mic's variability into con- 
sideration. The star regularly fiares in the X-ra y and EUV regions of t he spectrum (for 
a review with application to grain dynamics, see lAugereau fc Beustll2006l ). Periodic varia- 
tions in the visible regions of the spectrum, thought to be caused by spots (Amy ~ 0.35 mag, 
P = 4.865 d: [Torres &: Ferraz Melldll973l : ICutispoto et al.ll2003l ). are also relevant. However, 
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as noted by iMOSl . for AU Mic these are likely not problematic since the contrast between 



spots and the photosphere is lower in the near IR. 

It is convenient to measure surface brightness in rectangular photometric apertures, as 
are used in the calculations detailed in §3.21 Because of the difficulty in using coronagraphic 
observations to infer the off-axis (unobscured) PSF, we apply an aperture correction based 
on the enclosed energy in this aperture by using the unobscured observations of photometric 
standard stars as PSF references. In each band, we average the enclosed energy over az- 
imuthal rotation of the reference PSF to simulate the final roll-subtracted image's PSF. The 
aperture corrections we derived from these enclosed energy measurements are susceptible to 
variability in the AO PSF because the reference PSFs are measured non-contemporaneously 
from those of AU Mic. To estimate the systematic errors in the overall fiux levels, we exam- 
ined the random error in the enclosed energy of the reference PSF exposures in each band. 
The largest fluctuation on a single night was 16%. Although the combination of data from 
both nights will reduce this uncertainty (< 11%), the potential for unmeasured changes in 
the PSF from the reference star to science measurements remains. With caution in mind, we 
estimate the uncertainty in flux calibration error in the flnal near-IR images at 20% and note 
that, in general, unsensed fluctuations will tend to affect the shorter-wavelength observations 
to a greater degree. This systematic calibration uncertainty will affect the absolute levels of 
the surface brightness proflles we measure in §3.21 and depending on the degree in correla- 
tion between errors in calibration in different bands, the disk colors ( §3.3p . Measurements 
of color gradients will not be affected by this type of error. 



2.4. HST Imaging 

In order to compare disk images over a wider range of wavelengths, we reduced and 
analyzed data previously obtained with the Hubble Space Telescope Advanced Cam era for 



Surveys in the F606W flher (Ac=606nm, AA = 234 nm) on 2004 April 03 flKOSh . We 
acquired the flatfleld images of AU Mic and the PSF reference star HD 216149 from the HST 
OPUS pipehne. These data were further calibrated by dividing by an appropriate spotflat 
and multiplying by a pixel area map. A final image of AU Mic was constructed by averaging 
three frames of 750 seconds integration each, with appropriate filtering to reject cosmic rays. 
The final image of HD 216149 was constructed by averaging eigh t fra mes of 225 seconds 



each. We subtracted the PSF subtraction in a manner described by lK05l . There was no need 
to apply our specialized roll-subtraction technique (Appendix to these data because of 
the stability of the ACS PSF. The PSF-subtracted AU Mic image was then corrected for 
geometric distortion using a custom IDL routine (J. Krist, private communication). We 
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used the Tiny Tim PSF mo Krist &: Hookl |200J) to compute the enclosed energy in 
rectangular apertures to enable comparison of surface brightness profiles with our near-IR 
data. 

We used the same method for computing midplane surface brightness profiles described 
previously ( §2.31) . We checked the con sisten cy of our F606W midplane surface brightness 
measurements against those reported in lKOSl 's Fig. 4 (using their 0'.'25 x 0'.'25 apertures) and 
found that our measurements of the s ame data were uniformly 0.20 ±0.05 mag brighter. 
The discrepancy is due to the fact that iKOSi divided their images by a factor of 1.124 when 
correcting their pixel areas from 28 x 25 mas pixel"^ to 25 x 25 mas pixel" ^ (J. Krist, pri- 
vate communication). However, the photometric calibrations, both in the image headers 
and produced by SYNPHOT, assume that the final undistorted image is processed using 
DRIZ ZLE. The drizzled data product is similar to the manual calibration performed above 
and in lK05l . except that it does not include a uniform scaling of the i mage by a factor of 
1.124. To correct the erro r, the disk brightness values reported in lK05l should be increased 
by this factor. We adopted lK05r s value of 8.63±0.03 mag for the stellar brightness. Aperture 
size aside, our measurement methodology of the F606W midplane surface brightness profile 
presented in §3.21 produces results consistent with 'k05I after scaling their surface brightness 
by a factor of 1.124. The stability of the HST/ ACS PSF ensures that the flux calibration 
uncertainty, set by the 0.03 mag uncertainty in stellar brightness, is much lower than those 
of the near-IR bands. The stellar brightness in F606W was measured in unocculted images 
at the same epoch, so our measurements are unaffected by the star's variability. 



3. RESULTS 

3.1. Disk Morphology 

The PSF-subtracted images are shown in Figures [2] and [31 As observed by other authors, 
the disk has a near-edge-on morphology. The data in each band of Figs. [2] and [3] has been 
divided by the corresponding stellar brightness, such that the resulting colors trace the 
relative scattering efficiencies of the dust. We find the disk decreases in brightness relative 
to the star with increasing wavelength, up to i<''-band, indicating an overall blue color. We 
do not detect the disk with strong confidence in L' . 

The disk midplane shows vertical structure (Fig. |4]). We fit a Gaussian function to 
the vertical profile of the disk as a function of projected distance. We subtract this profile 



f 



http : //www . stsci . edu/ sof tware/tinytim/ 
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from the image, scale the residuals by a smooth fit to the amplitudes of the Gaussians, 
and display the results in the left panels of the Figure. This procedure clearly reveals the 
location of a sharp midplane. Broader features also show variation — the width of the 
best-fit Gaussians is not constant with projected distance from t he star. This effect has 
been noted and studied by other authors (IKOSl : iGraham et al.l 120071 ). Capturing the vertical 
structure of the disk is important for detailed physical models of the system. However, here 
and in ^ we restrict ourselves to considering only the most basic disk properties and defer 
detailed two-dimensional modeling for future work. We will discuss the disk substructure 
(brightness enhancements and deficits) in §3.4[ 



3.2. Surface Brightness Profiles 

Midplane surface brightness profiles are useful metrics of disk structure. These profiles 
average over the vertical extent of the observed scattered light image, which is naturally 
integrated along the line of sight. This facilitates comparison with scattered light models 
by reducing the brightness distribution to one dimension. We note that this convenient 
technique can fail to capture variation in the projected vertical extent of the disk, though we 
do follow the vertical centroid as a function of projected position. Prior to computing our 
profiles, we fit a spline to the vertical midplane position along each ansa. We used O'.'l x O'.'S 
photometric apertures centered on these positions in unsmoothed images to extract the 
photometry (Figure [5]). The width of this aperture is chosen to provide sufficient spatial 
resolution along the disk midplane, while the height is sufficien t to capture several vertical 



FWHM of the inner disk. As reported by iGraham et al.l (120071 ). the apparent F606W disk 



thickness increases outward, reaching the O'.'S aperture height at ~60 AU. Beyond this point, 
these photometric apertures do not perform a vertical average; rather, they sample the 
midplane brightness. 

Our PSF-subtracted L' data exhibit very low S/N structures along the disk plane in L' 
at separations of ~ 1" from the star, though these may result from PSF subtraction errors. 
In Figure [5l we show upper limits on the disk brightness in this band. 

The uncertainties in Fig. [5] represent the random measurement errors. At each radius, 
the contribution from residual speckle and background noise sources was estimated from the 
standard deviation of photometry in apertures placed in an annulus, excluding the locations 
near the disk. The photon noise from disk photometry is also factored into the random 
errors. The inner and outer edges of the profiles are set by the requirement that S/N > 1. 
To estimate the contribution of PSF subtraction to the systematic errors, at each radius we 
compared the mean values of the off-disk apertures (placed in an annulus) to the random 
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error estimated from their standard deviation. These are generally less than one, suggesting 
that any bias in the profiles is smaller than the estimated random errors. Another possible 
mechanism for systematic error arises in the near-IR data because the PSF estimates are 
derived from the target exposures ( §2.21 and Appx. If there is insufficient field rotation 
in the images, light from the inner disk may be subsumed in the PSF estimate, resulting in a 
self-subtraction of the disk. Nearly all of our datasets have APA > 30° (Table [1]), sufficient 
to exclude this as a source of significant error. However, the 2004 Aug 29 J-band data have 
APA ~ 17°, which corresponds to one disk midplane FWHM at ~ 8.5 AU. The 2004 Aug 
30 data have much more field rotation, so the effect on the inner region of the composite 
J-band surface brightness profiles is likely to be small. Finally, the uncertainties in Fig. [5] do 
not include potential errors due to fiux calibration, which may be as large as 0.2 mag for the 
JHK' bands ( §2.31) . We conclude that, relative to the random errors, the systematic errors 
are unlikely to affect the shape of the surface brightness profiles at separations greater than 
8.5 AU. 



By comparing our near-lR data to those of |L04J and |M05| . we confirm a break in the 



overall midplane surface brig htnes s profile around 30-35 AU. A break at ~ 15 AU in the 



F606W data was reported by lK05l . though the innermost points in our near-IR profiles are 



consistent with both a slight fiattening and no break at all. 

Broken power laws provide a compact description of the observed surface brightness 
profiles. Using fiux-based units rather than magnitudes, we fit a function f{b) oc to the 
midplane surface brightness profiles in Fig. [5l and report the resulting indices a in Table O 
We scale the formal l-cx random errors by \/xl when this quantity is > 1 to partially account 
for the ill-fittedness of a strict power-law to profiles which exhibit substructure. Errors in 
subtraction may be correlated because some PSF structures have significant radial extent. 
We ignore measurement covariance in the analysis, and thus potentially underestimate and 
the quoted errors. We chose to fit over two domains in projected separation: an inner region 
of 15.0-32 AU, and an outer region of 32-60.0 AU. The innermost boundary was chosen to 
mitigate possible biases from systematic errors. The outermost boundary is an upper limit, 
and is more precisely set by the availability of data in Fig. [51 The break between the two 
regions was chosen to correspond to the kink in the midplane surface brightness profile of 
the SE ansa, corresponding to feat ure C ( §3.4p . The power-law fits of our i/-band data are 



somewhat consistent with those of [LOJ, who measures for the NW (SE) ansa a = 1.4 ± 0.3 
(1.0 ±0.3) over 20-35 AU and 4.4 ±0.3 (4.4 ±0.4) over 35-60 AU. Using the same method as 
used in Table[2]on these ranges (but different aperture sizes and positions), we obtain 1.6±0.2 
(1.5 ± 0.1) and 4.4 ± 0.2 (3.4 ± 0.2), respectively. The values we compute for the SE slopes 
are stee per in the inner region, but shallower outside. Our measurements are consistent with 



those of M05 over similar domains. Those authors measure a = 1.2 ±0.3 over 17-33 AU and 
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4.0 ±0.6 over 33-60 AU. Using these ranges and averaging over ansae (but different aperture 
parameters) we obtain 1.4 ±0.2 and 3 .8 ± .2. We have significant differences in our F606W 
power-law indices relative to those of iKOSl . The authors, using apertures that are 0'.'25 tall 
in the direction perpendicular to the midplane, calculate a = 1.8 over the domain 15-43 AU, 
and 4.7 in the region beyond 43 AU, both of which are somewhat steeper than our findings 
of 1.6 ± 0.1 and 4.1 ± 0.2 (measuring from 15-43 AU and 43-70 AU with our apertures). 
However, measuring with O'.'l x 0'.'25 apertures over the region 43 -60 A U and averaging the 



fit indices over both ansae, we find a = 4.7 — consistent with K05. The differences in 



these comparisons serve to underscore the sensitivity of such power-law measurements to 
methodology. These differences may arise from PSF subtraction residuals, aperture sizes, 
aperture center locations, and the fit domain. Inferences between different datasets should 
be based on consistent power-law fits. 



3.3. Disk Color Variation 



The blue color, relative to the star, of the scattered light was first reported by IK05 



It becomes increasingly blue at larger radii {F435W-F814W = 0.2-0.5mag from 30-60 AU). 



The disk also scatters blue when comparing visible data to the i/-band measurements of lL04 



and lM05l . The blue color of the scattered light extends to our JHK' observations. Figure [H] 
shows the color as a function of projected radial position along the disk midplane. The stellar 
contribution to the apparent disk color has been removed in order to highlight intrinsic grain 
scattering processes. Systematic uncertainties in calibration are not included, though as 
discussed in §2.3H2.4l these are expected to be < 0.3 mag for K'-H and J-H, and < 0.2 mag 
for F606W-H. In the inner disk (< 35 AU), the observed color gradient is consistent with a 
fiat profile, while the outer disk is increasingly blue with stellocentric distance. 



3.4. Disk Substructure 



The AU Mic debri s disk is known to exhibit non-uniformity in its midplane surface 
brightness distribution (L04; Kod : Im05| ) . In order to remove global structure and highlight 



localized variations, we divided each brightness distribution by a smooth function in each 
band. We obtained weighted fits of spline functions to the NW and SE ansae's midplane 
surface brightness profiles ( §3.2p and averaged the results. The order and smoothness of 
the splines used to enhance the data have direct bearing on the spatial frequency content 
of the resulting map and may affect the observed positions and brightness of disk features. 
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We used the curf it routine of FITPACKg to fit a cubic spline to the data in Fig. [5l with 
smoothness s = m + ^/2m, where m is the number of data points in the profile. Figure [7] 
shows an F606W image and a JHK' color composite that have been processed in this 
manner. The annotations give the locations of substructural features (A-E). To allow for 
more quantitative comparison of disk features, in Figure [8] we show surface brightness profiles 
processed in the same manner. In order to highlight the structure of midplane brightness, 
these profiles are computed with O'.'l x O'.'l apertures. 

The imperfect removal of time-variable features in the PSF may introduce systematic 
errors which can masquerade as disk structures , un d erscor ing the value of independent ob- 



servations when identifying particular features. iLiul (120041) ide ntified several features of the 
disk substructure, which were also observed by iKOSi and iMOSl . In independent observations 
such as these, any changes seen in the structures may be due to PSF-subtraction artifacts. 
However, given that AU Mic exhibits brightness variability due to starspots, it is also pos- 
sible that these spots differentially illuminate the disk such that the spatial distribution of 
disk brightness is also time-variabile. We also note that differences in feature position and 
brightness from different authors may arise via analysis methodology. As we have done in 
Figs. [7] & [HI it is desirable to highlight features by dividi ng th e data by a function f{b), 
where b is the p rojected stellar separation. For example, lL04j u sed f{b) = \b\~^ (see his 
Fig. 3) and iMQa used f{b) = b~^ (their Fig. 3b). In contrast,!^ use a fourth-order polyno- 



mial fit to the surface brightness profile for /(&), while we use a spline fit (detailed above). 
These more complicated functions can better capture the overall brightness profile of the 
disk, which does not follow a single power law. On the other hand, they can suppress the 
appearance of the broadest features (e.g. feature D in the NW ansa). 

The application of the spline fit to the F606W- and J H K'-handi. data allows us to 
eliminate analysis methodology as a potential source of variation in feature brightness and 
position. While there are minor differences in the apparent characteristics of features between 
the near-IR bands, distinct differences in feature location arise when comparing those seen 
in visible to the near-IR. 

Feature A, which is a brightness enhancement seen in both ansae, is ~ 1 AU further 
from the star in the visible (26 AU) compared to the near-IR (25 AU). There is also evidence 
for a change in the position of this feature in the near-IR; in particular, the J-band centroid 
is further out than in H or K' (cf. Fig. [8]). Feature B is a brightness deficit in the SE ansa 
that also exhibits a similar inward-moving centroid with increasing wavelength. Feature C is 
a SE brightness enhancement, broader than A, that shows a different position depending on 



See http : //www . netlib . org/dierckx7| 
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the waveband. In J and the feature centroid is at ~ 32 AU, while it peaks closer to 33 AU 
in F606W and K' . Feature D (37 AU) corresponds to a broad enhancement on the NW ansa, 
a narrower dip in the SE ansa brightness, and the location of the vertical displacement of 
the NW midplane. The posi tion of this feature is not significantly wavelength-dependent in 
these data. We confirm iKOSl 's identification of a dip in midplane brightness inward of feature 
A, which we deno te 'E .' The locations of the other features seen in F606W agree with the 
positions given by iKOSl . The i f -band positions and characteristics of features A-D compare 
favorably to the results of lL04 Together these stren gthen the evidence for ~ 1 AU differences 
in the positions of features A-C. In contrast, iMOSi measures positions for A and B slightly 
closer to the star (by 2-3 AU); however those data are of lower S/N and the differences may 
not be significant. A summary of these results is given in Table [31 

There are additional deviations in brightness th at are detected in multiple bands. The 
broad enhancement at ~ 46 AU in the SE noted by lK05l is also detected in all the bands 
we consider. The feature seen in the F606W data, characterized by a peak enhancement at 
15 AU and sharper drop at 12 AU, has corresponding characteristics in the near-IR bands. 
In JHK', the peak is seen at 12-13 AU and the inner cutoff at 11-12 AU. The outer extent 
of this feature decreases gradually from ~ 20 AU in F606W to 17 AU in K' . The peak 
enhancement in the SE at ~12AU in F606W also has corresponding enhancements at 12- 
13 AU in J and if; an examination of the if'-band surface brightness profile (Fig. [5]) suggests 
the enhancement may also be present at ~ 13 AU in this band, but is not visible in Figs. [7] 
& [S] due to the spline fitting process. 



3.5. Point Source Detection Sensitivity 



We do not detect any point-like sources in the disk midplane. We have developed a 
methodology for utilizing the artificial insertion of point sources into an image to measure 
the detection sensitivities in both the background and the disk midplane. We detail our 
technique in Appendix [Bl and present our sensiti vity limits in F i gure P We indicate the 
predicted brightnesses of model giant planets from Burrows et al. ( 1997I ). placing limits on 
the presence of young, massive companions. We note the caveat that at y oung ages, the 



emission from Jupiter-m ass planets may be sensitive to initial conditions (e.g. iFortney et al. 
2OO5I : iMarlev et al.ll2007h . 
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ANALYSIS 



Models of the grain size and space distributions that reproduce the scattered light and 
thermal emission probe the dynamical processes affecting t he dust distribut ion. Previous 
models of the AU Mic di sk ha ve attempted to inodel t h e SEP (Chen et al.l2005 ) and th e scat- 
tered hght profiles (e.g. 



KOSl). iMetchev et all ((20051), lAugereau fc BeustI (120061 and ISC06 

i 



tackled the task of simultaneously fitting both. Using contin uous p ower-law descriptions for 

I find they cannot simulta- 



MO^ 



the grains' space and size distributions (sub-/im to mm sizes 
neously account for the scattered-light color, the SEP, and the break in th e surface brightness 
profiles at ~ 35 AU. As an alternative to power-law density distributions, lAugereau fc Beust 



(120061 ) obtained surface density models via scattered-light profile inversion. They tuned the 
size distribution to simultaneously fit the scattered light and SEP. A distribution of silicate 
grains successfully matched the F606W-ha.iad profiles and SEP, and in a separate fit, the 
i/-band profiles and SEP could be matched using a minimum grain size that is 10 times 
larger than the visible case. This discrepancy implies that size distributions which match the 
scattered light colors underestimate the am ount o f large grains responsible for the sub-mm 
emission. A different approach was taken by lSC06l . who modeled the dynamical structure of 
the debris disk to match the F606W scattered light and SEP. For sufficiently large values 
of the stellar mass loss rate M, this model produces a blue color, with an outward blue 
gradient. 

These previous efforts highlight two essential results. First, the blue color of the scat- 
tered light suggests small dust grains scattering in the Rayleigh regime, i.e., submicron parti- 
cles. Seceond, the shallow slope of the long- wavelength SEP is best fit with lar ge grains of up 



to mm size (grains larger than this are weakly constrained). More recently, iGraham et al. 
(120071 ) demonstrated that the linear polarization of the F606W scattered light requires 
porous grains. These data are important because the degeneracy between scattering asym- 
metry and spatial distribution can weaken inferences of debris disk structure based solely 
on measurements of total intensity. As the structure, size, and compositions of the grains 
determine the scattering properties, the assumption of compact grain types in previous work 
warrants reexamination. 

Our goal is to find the simplest description of the grains and their distributions that is 
compatible with the variety of available observational data. In §4.1l we describe the methods 
and structure of our models, the process used to fit the models to the observed scattered 
light and SEP, and the properties of the resulti ng best -fit models. We then examine these 
results in the context of the dynamical model of ISC06l for a steady-state grain distribution 
produced by a ring of parent bodies ( §4.2p . 
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4.1. Dust Modeling 

4-1.1. Monte Carlo Radiative Transfer 

We model the debris disk around AU Mic using MCFOST, a Monte Carlo radiative 
transfer code that uses the Stokes formalism to treat the interaction between dust grains and 
photons to p roduce SED s , scat tered light images, and polarization maps. MCFOST is fully 



described in iPinte et al.l (120061 ). and we summarize its main features here. Monochromatic 
photon packets are emitted by the star and propagated through the disk. The optical depth 
through which the photon travels before it interacts with a dust grain is randomly chosen 
from an e""^ probability distribution. Scattered light and the thermal SED are computed 
separately. When monochromatic scattered light images and polarization maps are produced, 
photons are only allowed to scatter off dust grains, with a loss of intensity corresponding to 
the absorption cross-section of the grain. When the SED is computed, photons are either 
scattered or absorbed, depending on the local albedo; in the latter case, they are immediately 
re-emitted at a longer wavelength selected on the basis of the local dust temperature. Once 
photons exit the computing volume, they are stored in "reception captors," corresponding 
to specific inclinations. Maps at all inclinations are simultaneously created, but in this study 
we focused on the most edge-on captor, which includes inclinations ranging from 88?9 to 90°, 
believed to be appropriate for the inner disk of AU Mic. 

To calculate the thermal equilibrium of dust grains, and therefore SEDs, the dust prop- 
erties must be known throughout the electromagnetic spectrum, from the ultraviolet to the 
millimeter regime. It is not possible to describe the dust grains with a simple parametriza- 
tion based on the albedo and the phase function asymmetry factor, for instance, unless these 
are known at all wavelengths. Rather, MCFOST relies on Mie theory (i.e., the grains are 
assumed to be spherical or essentially randomly oriented), so the knowledge of the dust 
optical indices (the complex index of refraction) at all wavelengths is sufficient. Grains 
with complicated structure, such as porous aggregates, are approximated by spheres with 
an effective optical index at each wavelength. Mie theory is used to calculate the scatter- 
ing and absorption properties for these effective-medium spheres. In the case of AU Mic, 
small Rayleigh-scattering grains are indicated by the blue color and high polarization frac- 



tion (IGraham et al.l 120071 ). For composite grains, an effective medium approximation (de- 



pending on the constituent materials and me thod for computing the effecti ve indices) can 



be reasonably accurate in the Rayleigh limit (IVoshchinnikov fc Mathislll999l ). 



MCFOST was first developed to model gas-rich, optically thick disks surrounding 
T Tauri stars. Several features have been added to efficiently and correctly treat the case of 
optically thin debris disks like AU Mic. The first modification was to increase the computa- 
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tional efficiency of treating optically thin material. Rather than expend effort computing the 
fates of all random photons, the majority of which will not scatter off grains in the optically 
thin disk, we enforce the first scattering of each photon packet to occur within the disk. The 
expected non-interacting photons are not randomly generated; instead we analytically ac- 
count for the corresponding transmitted starlight. This ensures that all randomly generated 
stellar photons scatter at least once in the disk without energy loss. Second, in the absence 
of gas and given the low density of dust particles, the dust grains cannot be considered to be 
thermally equilibrated with one another. Rather, we compute a size-dependent temperature 
for the grains, each being in equilibrium with the surrounding radiation field. In the case of 
the AU Mic disk, this results in a significantly different shape of the SED in the mid-infrared 
due to large difference in temperature between grains of different size. Future high-resolution 
spectra in the 10 /xm region m ay measure features which can constrain the composition and 



size distribution of grains (e.g. iLi &: Greenberg||l998l : IChen et al.ll2006l ). 



4.1.2. Model Construction 



Disk models frequently assume power-law descriptions for both the geometry of the disk 
and the grain size distribution. This is unlikely to reproduce all observed aspects of the disk, 
such as the complex small-scale structure seen in the scattered hght ( §3.ip . It is nonetheless 
a valuable approach for constraining some of the main system parameters and testing simple 
hypotheses. Here, we attempt to reproduce the observations of the disk (SED and scattered 
light proffies) with a two-zone disk description. In a narrow inner zone, large grains would be 
present and account for most of the long wavelength thermal emission, whereas a much more 
extended outer region would contain small grains and be resp onsible for the scattered light. 



This i s a model that qualitatively matches models proposed by lSCOa and lAugereau fc Beust 



(I2OO6I ). where our inner annulus would represent the observable population of parent bodies 
and the outer zone the populations of small grains created by collisions of large bodies and 
swept out by pressure forces. 

In both regions, we assume that the surface density and the grain size distributions 
follow S(r) oc and dN{a) oc a'^da. The outer radius is fixed at 300 AU , as scattered light 
has been detected out to a sensitivity limit of 210 AU (IKalas et al.l 120041 ). The inner radius 
and transition radius between the inner and outer regions are free parameters. In addition, 
the surface density index, minimum and maximum grain size and distribution index, and 
the total dust mass in each of the two zones are free. We initially use a 7 = —7/2 index 
for the size distribution, which is c ommonly assuni ed for debris systems and is suitable for 
a steady-state coUisional cascade (IDohnanyil Il969l ). While both theory and observations 
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may suggest specific functional forms for the vertical density profile of the disk, we cede 
to our preference for simplicity and assume that the vertical density distribution follows 
a Gaussian with a fixed width of a = 0.8 AU and a flat radial dependence. In this case, 
the photometric apertures capture the majority of the scattered light flux during brightness 
profile calculation. The particular choice of vertical density function is unimportant as long 
as the scale is much less than the photometric aperture height ( §3.21) . We defer for future 
modeling the exploration of relationships between the form of the vertical density profile, the 
dependence of scale height with radius, grain scattering phase functions, and the geometry 
of the disk. 



The data for the SEP were firs t compiled by iMOSi and supplemented with new Spitzer 
measurements bv lChen et al.l (|2005l). We use a Tpff = 3600 K, \ogg = 4.5 NextGen model 
for the stellar photosphere (IHauschildt et al.lll999l ). using radius R = 0.88 i?© to match the 
observed stellar flux. For scattered light, we use the J-, H-, and K'-hand profiles from data 



presented here, as well as the F606W-ha.nd profiles from data initially presented by IKOSI and 
re-analyzed for this work in §3.2[ 



4.1.3. Fitting Procedure 



Thorough exploration of the available parameter space is a time-consuming task when 
producing a SED and multiple scattered light images. We therefore explored it by hand, 
iteratively narrowing down the possible values for each parameter. Our objective is not to 
find the best possible model, but to determine if there is at least one solution based on this 
simple two-zone model that can match th e data relatively well, whic h would in turn provide 
support to theoretical birth ring models (ISCOa : lAugereau fc BeustI 12006 ). We compare our 
results to these models in §4.2[ To further reduce computation, we have elected to fit the 
F606W- and i/-band profiles and the SED, and to later check for consistency with J- and 
i^'-band profiles and F606W fractional polarization curves of I Graham et al.l (120071 ). We 
calculate the model profiles with the same rectangular aperture sizes as the observed profiles 
in 3321(0'.' 1 X 0'.'5). 

We used a grain model that intrinsically includes a high porosity, a characteristic that 
provides b oth strong forward sc attering and polarization — nec essary features as demon- 
strated bv iGraham et all J2007h . The dust model developed bv iMathis fc WhiffenI Jl989l ) 
consists of a mixture of silicate, carbonaceous, and icy small elements combined into porous 
aggregates. With a^in = 5nm, Omax = 0.9 /xm, a size distribution index of 7 = —3.9, and 
a vacuum fraction of 80%, this model reproduces the interstellar extinction law from the 
ultraviolet to the near-infrared. We use the same optical constants but allow the size distri- 
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bution to vary. We fixed the inner zone's size-distribution index to 7 = —3.5. We searched 
for the power law distributions that would adequately produce the blue color (in terms of 
albedo x opacity) and roughly wavelength-independent phase function, so as to maintain 
roughly parallel surface brightness profiles from F606W to H. It was then possible to adjust 
the model parameters until a match to the SED and the scattered light profiles was obtained. 
We tuned the boundaries, grain size distributions, surface density indices, and total mass 
of each zone. For c omparison, we also fit a model using the optical properties of compact 
silicate grains (from lDraine fc Lill200ll ) using the same procedure. 



4-1 .4- Fit Results 

The parameters of our best-fit models are listed in Tables H] and [5l and the results are 
shown in Figures [10] and [111 Considering the F606W- and i/-band profiles and the SED, we 
obtain reasonable agreement with the predictions our two-zone model using porous grains: 
a region of large grains between 35-40 AU, and a region of smaller grains outward of 40 AU. 
These are slightly further out than the best-fit regions using compact grains, which are more 
isotropically scattering. We find that, if the particles in the inner region are porous, they 
must be in the mm-regime (sizes of 3-6 mm). Allowing the size distribution to encompass 
smaller grains results in an overly steep long-wavelength SED as well as a flux deficiency 
in the outer small grain region. Much larger bodies could be present yet undetectable; 
they would emit very little, even at mm wavelengths. For porous dust in the outer region, 
we require a distribution encompassing small grains (0.05-3.0 /im, dn oc a~^^^da), with a 
power-law index that is steeper than the collisional steady-state value of -3.5. 

The masses in the inner and outer regions of the porous grain model are around 1.0 x 
1O~^M0 and 2.3 x lO'^^M®, respectively. The surface density profile falls rapidly outside 
of the transition radius (index -2.5 to -3.0), as expected from swept up material; shallower 
indices yield an improved SED but insufficient scattered flux inside 40 AU and too-shallow 
profiles overall. Steeper indices yield an excessively precipitous F606W-hBxid brightness 
profile and excess flux around 50-100 /xm. In the inner region, the surface density can be fiat 
or increase with radius, peaking at the transition radius. This is only loosely constrained: 
the big grains never get very hot in this narrow annulus and there is little effect on the SED. 

We performed some exploration of the parameters to estimate uncertainties. The masses 
are relatively well constrained, to 50% or so, depending on the grain size distribution pa- 
rameters. Again, more mass could be hidden in the form of very large bodies in the inner 
region; the mass derived here is just what is needed to reproduce the long wavelength end 
of the SED. In computing the total intensity of scattered light, there is an anticorrelation 
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between amin and the region boundary at ~ 40 AU. We can adjust the minimum grain size 
by a factor of 2, and the ring radius by ±5AU and still achieve satisfactory F606W- and 
i/-band profiles and SED. 



As demonstrated by iGraham et al.l (120071 ). the fractional linear polarization profiles 



provide a complementary constraint on the optical properties of the small grains. The 
measured fractional polarization gradually rises to a plateau of ~35% at around 50 AU. Our 
best-fit models using compact and porous grains are comparable in fit quahty for the total 
scattered light intensity profiles and the SED; however the polarization curves in Fig. [TT] 
clearly rule out compact grains, which reach maximum fractional polarizations < 10%. The 
porous grains generally follow the measured trend, and reach a plateau of 30-35% by 40 AU, 
the starting location of the small outer grain region. 

This model has some shortcomings, in that it predicts excess flux around 60-70 /xm and 
a slightly too-shallow scattered-light profile in the H band outside of the transition radius. 
This latter point results in a predicted disk color that does not change outside of 1", as 
opposed to the observations. We also note that the model profiles show less agreement with 
the J- and if '-band profiles (Fig. [TTi) . The overall flux levels differ (though this may coincide 
with photometric calibration uncertainty; cf. §2.3p . and we overpredict the midplane surface 
brightness in the outer regions. For this analysis, we have avoided goodness-of-fit metrics like 
because our model, by design, will fail to reproduce substructure in the scattered-light 
profiles. Further, PSF subtraction residuals introduce correlations between errors in the 
measured profiles, which we ignore ( §3.21) . Nevertheless, we can get a sense of model fidelity 
by considering the long wavelength end of the SED (where the dust contributes). For the 
nine SED points in Fig. [10] at A > 10 /im, the best-fit compact grain model gives x = 8-6 and 
the porous model gives x = 10.9. The bulk of the deviation arises from the 60-70 /im region. 
A formal x^ analysis is not applicable, as we have also tuned our parameters to match the 
scattered light. In contrast to the porous grain model, the compact grain model underpredicts 
the 60-70 yum flux. With their higher average albedo, the compact grains absorb less energy 
than the porous ones (Table [5]), suggesting a better fit may be obtained with grain albedos in 
the 0.6-0.7 range. Still, the fit is satisfying, considering the simplicity of the model ( §4.1.2p . 
Presumably, considering non-power law prescriptions, other dust composition, non-spherical 
grains, and/or overlap between small and large grains would result in a better fit. One 
attractive possibility, that we do not explore here, is that ever-smaller grains are present as 
we move outside of the belt of parent bodies ( §4.21) . This could account for the increasingly 
blue color in the outer region as well as for the seemingly different transition radii obtained 
from the F606W- and H-h&nd profiles (~40AU vs. ~30AU). 
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4.2. Birth Ring Examination 



In the Solar System, the Asteroid Beh and Classical Kuiper Belt constitute rings of 
solid bodies that, through mutual collisions, act as sites of continual dust production. Several 
forces act on the newly liberated grains, and the resulting trajectories shape the overall struc- 
ture of the Sun's debris disk. A similar belt of parent bodies ha s been invoked to explain the 



break in the surface brightn ess profile of the (3 Pic debris disk ( iLecavelier Pes Etangs et al. 



19961 : lAugereau et al.ll2001l). Additional ri r igs have been prop osed to explain features seen 
in recent imaging fe.g. IWahhai et al.ll2003l: iTelesco et al.ll2005l ) and spatially resolved spec- 



troscopy (jWeinberger et al.ll2003l : lOkamoto et al.ll2004j ) of this system. An attractive feature 
of such models is that their dust distributions are in steady-state — we need not appeal to a 
rare (but recent) catastrophic collision between large (e.g. km-sized) bodies, such that fresh 
grains have had insufficient time to completely diffuse through the system. 

The breaks in the slopes of AU Mic's surface brightness profiles naturally raise the 
prospect of a ring of parent bo dies analogo us to the Classical Kuiper Bel t. The existence 



of such a belt was proposed by IMOSI . ISCOd and lAugereau fc Beust 



(120061 ). Here, we focus 



our attention on a detailed theoretical model, developed by SCod . which yields the steady 



state spatial and size distributions of dust grains in the disk. A critical parameter is the 
stellar mass- loss rate, M^, which governs the corpuscular forces on the grains. In addition to 
modifying the disk structure, the stellar wind is cru cial in determining the lifetimes of debris 
disks around late-type stars (jPlavchan et al.ll2005l ). Assuming the scenario of steady-state 
dust production by a ring of parent bodies is applicable, one may use a dynamical model to 
infer from observations. W e note that the theory was applied to the AU Mic disk prior 
to the observational findings of iGraham et al.l (120071 ). who showed that the disk's scattered 
light is strongly polarized, and that porous grains are required to fit the polarization profile. 
In this paper, we examine the disk structure predicted by the birth ring scenario in light of 
the currently available observational data. After a brief review of the predictions for compact 
grains, we probe modifications to the theory given the porous grains of our best-fit model 
and check the consistency of the theory's features with those of our Monte Carlo radiative 
transfer model. In this context, we appraise the predictive power of the observations to 
determine M^: via inferences of the size and space distributions of grains. 



4-2.1. The Theoretical Scenario 



We first recapitulate the basic physical arguments of the model put forth by lSC06l . Col- 
lisions occur within the ring of parent bodies, producing dust grains. The radial components 
of forces arising from stellar radiation and wind result in a pressure acting to push grains 
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away from the star, causing blow-out for sufficiently small grains. In contrast, the tangential 
components of these forces act as a drag mechanism, which for radiative forces is the well- 
known Poynting-Robertson (PR) drag. Grains dominated by drag forces spiral in toward 
the star. The competition between these forces and the relative frequency of destructive 
grain c ollision s determine the dynamical structure of the disk. In disks dominated by drag 



forces, ( iSCOGP s "type A"), the action of PR and corpuscular drags fill the region interior to 
the birth ring with grains that avoid destructive collisions. In disks dominated by collisions 
("type B"), this interior is essentially empty, as the pulverized remains of colliding grains 
are quickly swept outward. In both cases, the region exterior to the birth ring is largely 
populated by the small, tenuously bound grains which follow elliptical orbits with periastra 
near the birth ring. In type B disks, these grains are just larger than the blow-out size of 
stellar wind and radiation pressure. In type A disks, the smallest grains are pulled inward 
by corpuscular and PR drag, and the peak of the grain size distribution in the outer disk 
corresponds to a size significantly larger than that of blow-out. 



SC06I developed analytical and numerical descriptions of the results of these processes 
on the size-space distributions of dust grains, and then applied this theoretical framework to 
the case of AU Mic. They postulate a belt of parent bodies at 43 AU, a location determined 
by the surface brightness profile break in HST imaging. By sim ultane ously modeling the 



scattered light profiles {F606W -hand from iKOSi and if-band from iMOSi ) and the SED with 



compact spheres of pure water ice, they determine that type B conditions hold and place 
limits on M^.. With these assumptions on grain type, there are two pieces of evidence for an 
inner disk devoid of in-spiraling grains, which is expected for type B conditions: (1) a lack of 
photospheric excess in the 10 /im region of the SED, and (2) the shallow slope of the surface 



bright ness profile at projected distances interior to the proposed ring. lAugereau fc Beust 



(120061 ) argue for a similar dynamical scenario for explaining the observed structure, where 
corpuscular pressure forces sweep grains outward. Considering the pressure from the stellar 
wind (and additional radiative enhancement by stellar fiares), they determine that radial 
forces can sufficiently diffuse small grains into the outer regions to account for the observed 
profile break, in analogy with the dominance of radiation pressure in (3 Pic's disk. 



4-2.2. Potential Ramifications of Porous Grains 



It is important to examine the birth ring scenario in view of t he new constraints on scat- 
tering properties from recent imaging of the disk in polarized light (iGraham et al.ll2007l ). The 
polarization data strongly indicate the presence of porous grains scattering in the Rayleigh 
regime, reaffirmed in our modeling ( §4.1.41) . Such grains are largely forward-scattering and 
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produce the required high peak polarization fraction. Since scattering asymmetry directly 
affects the inferred spatial distribution of grains, previous models' use of spheres of solid 
material must be re-evaluated with this new evidence. We also stress that, as noted in §3.2^ 
analyses involving surface brightness profiles require consistent methodology across all data 
set s. Fol lowing these points, it is appropriate to both check the consistency of the model 



by ISCOGI with our porous grain Monte Carlo model, and to evaluate the strength of such 



models for inference of AU Mic's stellar mass-loss rate. 

In approaching this problem, it is apt to contrast the relevant properties of the com- 
p act and porous gram type s we use in the models of §4.1[ As expected, the porous grains 



of iMathis &: Whiffen (|l989f) are less de nse than the compact grains, at p = 0.5 gem (com- 



pared to 2gcm~ for lDraine fc Lill200ll grains). At a given size, the porous grains also couple 
less efficiently to the radiation field than their compact counterparts. Taking a 0.05 /xm grain 
as an example, the effective cross section for radiation pressure (Qprvra^) is more than an 
order of magnitude smaller for a porous grain: Q^^ ~ 2 — » 0.05 when changing from compact 
to porous. This coupling inefficiency is similarly manifested in a decrease in the effective 
cross section to scattering, QscaTra^, for porous grains. 

Debris disk structure is commonly parameterized by the vertical optical depth. In the 
following analysis, we will refer to both the vertical optical depth to scattering, rj^^, and the 
geometric vertical optical depth, . The former quantity is useful for describing the 1-d 
brightness profile, while the latter governs the coUisional timescale. For a given surface mass 
density S, the depth depends only on the grain density and size distribution dN/da, 
while T^'^ additionally depends on scattering efficiency: 

rfo^r) oc ^ / (na') ^da, (1) 
p J da 

y(r) r dN 
rr(r) [QsM^a'] —da. (2) 

For a given grain type and size distribution, Tf^° and can differ by orders of magnitude. 
As an example, consider the outer region of porous grains in our model from the previ- 
ous section. While the size distribution is steep {dN/da oc a~^'^), because Qsca(l A^ni) S> 
<5sca(0.05 pm), micron-sized grains contstitute the bulk of Tgca despite their relative scarcity. 
The geometric depth Tgeo, on the other hand, is dominated by the smallest grains, resulting 

111 ^geo ^sca- 
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4-2.3. Disk Structure Analysis: Theoretical Considerations 



The models of ISC06l delineate three regions of water ice grains: (1) a birth ring, popu- 
lated by parent bodies of sizes up to ~10 cm, (2) an outer region made up of smaller grains 
(c^min ^ 1 /^m) on loosely bound orbits, and in type A disks, (3) an inner region of grains, 
<//m size, spiraling inward as a result of corpuscular and PR drag forces. Grains smaller 
than the stellar wind and radiation pressure blow-out size, amin (~0.1;um), have negligible 
contribution to the disk. Our modeling of the dust distribution with porous grains ( §4.11) 
shows that a two- zone architecture (lacking an inner region of drag-dominated grains) is con- 
sistent with the observations and, qualitatively, with the architecture predicted by physically 
modeling dust generated by a belt of parent bodies. 

In the birth ring model, how does the expec t ed gr ain size distribution change in light 
of the differences between the iMathis &: WhiffenI ( 119891 ) grains and water ice spheres? We 
first consider abiow? which is determined by the combined action of radiation and stellar wind 
pressures. The density decrease of porous grains relative to spheres tends to increase Obiow 
However, this effect is offset by the decrease in the effective cross section, parameterized 
by Qpr (a factor for 40 lower for a = 0.05 /im grains), resulting in a n et dec rease of porous 



grains' blow-out size relative to their compact counterparts. Rewriting [SCOg's Eqs. 7 and 8, 
the blow-out size assuming constant pressure is 



Oblc 



87r GM^cp 



wind" 



(3) 



We follow ISC06I in adopting t^^nd = 450 km s^^ and Mq = 2 x 10^^'^ Mq yr~^ We choose 
L = 0.12 and assume that the cross section to wind (Qwind^ra^) is unchanged (Qwind = 1) 
for the porous grains. Together these yield Obiow = {0.029, 0.047, 0.23, 2.1 } fim for = { 1, 
10, 10^, 10^ } Mq. The minimum grain size in the outer region of our disk model (0.05 /im) 
is larger than Obiow for the weaker stellar winds, which suggests <C lO^M©. 

Next, we consider the implications of our porous grain model on the amount of disk infill 
(type A vs. B) predicted by the birth ring theory — a process regulated by the collision 
rate, which depends on t]^° (Eq. [1]). In our model for the scattered light at tbr = 40 AU, 
the smallest grains contribute the bulk of the total geometric cross section. The abundance 
of these grains increases tJ^° at the birth ring relative to the ISCOGI value by a factor of 10, 
to Tj^°{rBR) = 4 X 10~^. This shortens the timescale for destructive collisions of the average 
particle, assuming the specific collisional energy is unchanged. The coUisio nal tim escale is 
balanced by the drag-induced infall time at a grain size defined by flbreak (cf. ISCOGI Eqs. 36- 
37). Grains near this size constitute the bulk of T|^^°(rBR), and in equilibrium, grains of sizes 
c^biow < a < ttbrcak AH the interior (r < tbr) . Like the case of compact ice spheres, for porous 
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grains we find that abreak is the same order as abiow: ctbreak = {1.005, 1.03, 1.08, 1.1 } abiow for 
a peak g eometric vertical optical depth at tbr and the same stellar wind values. In contrast 



to iSCOy's results for the compact grain case, we find that (abreak — awow) < awow for porous 
grains in all wind cases, including stronger values. The implication of this result is that no 
amount of corpuscular drag will produce a significant fraction of grains with orbits crossing 
the inner region, favoring a clear inner disk (type B). 



4-2. 4- Disk Structure Analysis: Observational Limits on Infill 

As a complementary approach to examining the birth ring scenario, we can, for the 
moment, set aside the final conclusion of the previous section and use our model from §4.11 
to test allowed small-grain infill. When using MCFOST to model the observed disk with 
porous grains, can we make a distinction between the conditions of type A and B? By 
definition, our two- zone model lacks grains inward of the putative birth ring, consistent with 
a collisionally dominated type B disk. To test whether the system could be consistent with 
type A, we now consider adding a third region of grains. If a significant mass of small 
grains are present inward of ~35AU (type A), then these warmer grains would increase 
dust emission in the 10-20 /im region of the SED and add to the scattered light of the 
system (Figs. [TUHTT]) . We note that, compared to more isotropically scattering spheres, the 
strong forward scattering of porous grains allows us to hide more mass in this region without 
strongly affecting the total intensity of scattered light. The smallest projected separations, 
where the added contribution from forward scattering would be maximal, are lost in the 
glare of the central star and are behind the occulting spot. 

The additional zone of porous grains has a flat S(r) oc r° density profile and r anges 



from 15-35 AU. The inner cutoff corresponds to the estimate of the ice boundary by |M05 



(13-15 AU). As a general feature of the birth ring conditions modeled by ISC06l . the outer 
region's grains have amin ~ c^break ~ c^biow; therefore we model a single-sized population of 
a = 0.05 fim grains in this new inner region. The middle and outer regions' dust distributions 
were fixed to the best-fit parameters (Table H]), and we adjusted the mass in the new inner 
region in ~ 0.5 dex steps until the scattered light profiles (including fractional polarization) or 
SED showed significant deviation. We found that the SED first showed noticeable deviation 
at M = lO-^M®; this is an upper limit to the mass of grains in this population allowed by 
the observations. 

We can express this limit as either the geometric or scattering vertical optical depths 
relative to the peak depth in the birth ring, ( = T±{r < tbr) / T±{rBR) ■ The geometric vertical 
optical depth is more closely related to the dynamical structure, while the scattering depth 
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is traced by our imaging observables. As noted in §4.2.21 these geometric and scattering 
optical depths can be quite different. We find r|'^°(r < tbr) < 2.5 x 10~^, giving ^gco < 0.6. 
While our limit to Cgeo is order unity, the strong decline in Qsca. with grain size results 
in a more stringent limit to Csca for 0.05 /xm grain filling. An inner mass of 3 x 10~^M^, 
easily r uled out by our analysis, corresponds to Csca ~ 1/400. This is the same order as the 



limit of iGraham et al.l (120071 ). who determined Csca ^ 1/300 inward of 40 AU using a similar 



method (though without SED modeling). 

While the limit for (^gca implies an inner disk relatively free of scattering grains, the rele- 
vant quantity for the dynamical structure of the disk is Cgeo- The weakness of our geometric 
optical depth ratio limit implies that we cannot yet determine type A vs. B dynamics or 
make inferences of M^, by matching the expected disk structure to current scattered light 
and SED data. Rather, perhaps the strongest direct observational evidence for a type A disk 
would come from resolved thermal emission from this zone. 



4-2.5. Future Considerations 



When considering more rigorous tests of birth ring theory, several enhancements can be 
made to our simple disk model. Faithful models must not rely on an artificial separation 
between the site of dust production (inner, large-grain region) and the small grains, as we 
have done in §4.1[ ISCOa also predict different asymptotic power-laws for the outer region's 
geometric vertical optical depth in type A and B disks (-2.5 and -1.5, respectively). However, 
the corresponding asymptotes in scattered light profiles are not reached until distances of a 
few 100 AU, ruling out reliable comparisons with current data. Somewhat more challenging 
aspects of the t heory are grain size distri butions that do not follow power-law description (cf. 
Fig. 3of|sC06|). The prediction of ISCOGi is an excess of grains above the Dohnanyi collisional 
cascade {dn oc a~'^^'^da) for sizes near auow Incorporating a population of grains at this size, 
in addition to the p owe r-law size distribut i on, is a possible avenue for future significance 
testing. Both ISCOGI and lAugereau &: BeustI (120061 ) predict diminishing minimum grain size 
with increasing stellar distance (due to ever-more tenuously bound grains as a result of stellar 
wind and radiation pressure), which can effectively model the observed blue color gradient 
in the outer region. Future scattered light modeling may incorporate a function amin(r) in 
the outer zone. 
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DISCUSSION & CONCLUSIONS 



Debris disks are a long-lived phase in the evolution of circumstellar material, and they 
provide indirect probes for the presence of planets and planetesimals. Resolved imaging is 
crucial for determining the structure and dynamical history of the dust. In our near-IR 
imaging, we measured the structure of the disk see n in scatt ered light . We confirmed the 
overall architecture of the disk seen in H band by |L04| and iMOSl . and presented the first 
images seen in J and K' ( §3.ip . We represented the global disk structure with a power-law 
description of the midplane surface brightness, and noted that discretion must be used when 
comparing between profiles measured in different manners ( §3.21) . Scattered light features 



on both large and small scales are largely consistent across the bands studied here. The 
colors we measured in the surface brightness profiles ( §3.31) allow inferences of the scattering 
properties, which are determined by the composition, sizes, and structures of dust grains. 

We used a model for the space and size distributions of the grains to reproduce these 
profiles and the SED ( §4.11) . The simple, two-zone description of power-law grain distribu- 
tions cannot faithfully reproduce the observations with compact silicate grains. Rather, 
porous grains of dirty ice , whose presence is indicated by recent polarization measure- 
ments ( iGraham et al.ll2007l ). have optical constants with which we can match such an archi- 
tecture to the data with reasonable fidelity. The key degen erac y is b etween the scattering 
phase function and the density distribution. Models of iKOSi and iMOSi place grains which ex- 
hibit moderate scattering asymmetry in the inner regions of the disk. The porous grains in 
the outer zone of our model, which are larger and more forward scattering than the compact 
grains used in previous efforts, fill in this region of scattered light. As such, our modeling 
does not require grains inward of 35 AU to reproduce the midplane surface brightness profiles. 
The inner zone (35-40 AU) of our model is made up of few-mm-sized grains, which reproduce 
the far-IR end of the SED. The outer zone grains of our model are smaller (0.05-3.0 /xm) 
and produce the majority of the mid-IR emission and the scattered light. 

The blue color of the s cattered light disk is relatively rare among resolved systems (cf. 
Table 1 of lMeyer et al. I l2007h . The succession of modeling efforts have shown that a significant 
population of submicron grains reproduces the disk's blue scattered-light color, and that the 
presence of such grains is compatible with interpretations of the disk's dynamical structure. 
As other authors have noted, stellar properties play a major role in the removal of submicron 
grains. As an M dwarf, AU Mic has a uniquely low radiation pressure relative to other 
stars with scattered-light debris disks. For sufficiently low values of M^:, we find awow is 
submicron and grains will scatter blue. While stellar properties can reasonably be used as 
a predictor for the presence of submicron grains in debris systems, grain properties can also 
affect the blow-out size. As we noted in §4.2.2[ the smallest of our porous grains have much 
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smaller Qpj- relative to their compact counterparts of the same size. The blow-out sizes we 
calculate for por ous gra ins are up to a factor of eight smaller than those of compact grains 
(as computed by ISC06l ). When blue scattered-light colors are dete cted around earlier-type 



stars with stronger radiation pressure (e.g. HD 32297 and HD 15115; lKalasll2005l : iKalas et al. 



20071). inferences of submicron grains may constrain their porosity and composition (see also 
Artymowica Il988l : iGrigorieva et al.l 120071 ) . 



The study of debris disk evolution hinges on accurate measurements of stellar ages and 
dust masses. The total observed mass of the disk in our model (1.0 x 10~^ M®) is dominated 
by the mm-sized grains in the inner zone. Note that this is a factor of few times smaller than 
the mass we derive for compact grains (2.7x 10~^ ^e)- Though t he equilibrium temperature s 
for the two grain types differ, this is primarily an opacity effect ( iVoshchinnikov et al.ll2006l ). 
The decrease in density increases the opacity, which in turn lowers the mass needed to 
match the sub-mm fluxes. Since these grains contribute little to the scattered light, we 
have no information about their porosity (and correspondingly whether any compaction has 
occurred). Therefore, in addition to the normal caveats about dust masses derived from sub- 
mm fluxes (e.g. insensitivity to larger bodies), derived debris disk masses are additionally 
uncertain by a factor of a few depending on the grain porosity. The same opacity caveat 
applies to the inferred sizes of the grains that dominate the sub-mm emission. In the inner 
region, our porous grains are a factor of three larger than the equivalent best-fit compact 
ones — the increased opacity means a larger size is needed to match a given temperature. 

Nevertheless, we can compare our dust masses to those of other authors. Our estimate 
is similar to the rn ass estimated by a single-temperature fit to the l ong- w avelength SED 
bv lLiu et al.l J2004h . who found 1.1 x lO^^M® of 40 K dust at 17 AU. Im05I found a similar 



mass as iLiu et al 



, for sing l e-zon e power-law distributions of size and density of ISM-like 



grams. lAugereau fc BeustI (120061 ) calculate a mass of 7 x lO^'^ by adjusting the grain 
size dist ributi on in their scattered light model in order to match the long-wavelength SED. 
Finally, ISCOGI argue for parent body sizes up to 10 cm, the maximum size participating in 
a steady-state collisional cascade. For their model's size distribution, the mass of parent 
bodies is ~10~^M^. However, current observations are not sensitive to grains larger than 
a few mm, so the maximum observed mass in this model is significantly lower. For siz e 

1 /2 

distribution index 7 = —7/2, the mass M oc Omax, and the observable portion of iSCOGr s 
mass is a few times smaller than that of our work. 

The interpretation of debris disk observations for inferring the presence of planets and 
the history of the grains rests on sound physical modeling of grain dynamics. The overall 
shape of the midplane surface brightness profile distribution is consistent with steady-state 
dynamics arising from a ring of parent bodies ( §4.2p . With models of grain dynamics, we 
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might constrain M by examining the density of small grains in the inner disk (< 35 AU). We 
found that for the porous grains in our model, an inner region would be populated primarily 
by grains just above the blow-out size for a wide range of M^. Further, the constraint of the 
presence of grains in an inner region is set by the mid-IR SED, and we can only say that 
the geometric vertical optical depth in this region is at most the same order as the peak 
depth in the birth ring, consistent with both drag- and collision-dominated (type A and B) 
disks. However, the presence of small (0.05 /xm) grains in the outer region implies that the 
mass-loss rate is low enough such that these grains survive blow-out, i.e. M^, ^ IO^Mq. 

The overall disk architecture inferred from these models has implications for the sub- 
structure seen in scattered light. In §3.41 we compared our measurements of features A-E 
with the results of other authors, giving further evidence that they are not artifacts of stellar 
PSF subtraction. Additional brightness enhancements located 11-15 AU from the star have 
also been identified, though additional observations of the disk at small angles are needed 
to improve confidence. The appearance of such substructural features is determined by the 
grain distribution; since the bulk of the scattered light comes from small grains beyond 
~ 40 AU, the features at projected distances inward of 40 AU must arise from azimuthal 
perturbations in the dust distribution beyond this distance. Some of these features appear 
to exhibit a wavelength-dependent position. (We note the position of feature D, which is 
closest to 40 AU, does not appear to do so.) While the overall position of a feature may be 
used to trace a perturbation in the density distribution in the outer disk, the wavelength 
dependence of the position may constrain additional spatial variation in the grain size dis- 
tribution. In a similar vein, any vertical structure associated with a given feature (Fig. H]) 
can be tied to the vertical structure of the outer disk. 

In summary, we: 

1. Demonstrate a new roll-subtraction technique that attempts to mitigate AO PSF vari- 
ability (§221 Appendix ED, 

2. Detect the AU Mic debris disk in Ji/ if'-band imaging and place an upper limit on the 
L'-band brightness distribution ( §3.11) . 

3. Confirm the blue color of the disk and measured a blue color gradient outside of the 
transition radius ( §3.2113.31) . 

4. Place detection limits on point sources in the disk midplane ( §3.51) and document 
a technique for determining point-source detection sensitivity in the disk midplane 
(Appendix [B]), 
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5. Verified the presence of substructure in tlie inner disk, and sliown tliat some features 
exhibit shght variation in their positions with wavelength ( §3.41) . 

6. Demonstrate the apphcabihty of a simple two-zone model which simultaneously fits 
the scattered light profiles and SED ( §4.ip . 

7. Find that compact silicate grains cannot be used in our simple model, while porous, icy 
aggregates of silicate and carbonaceous grains can reasonably account for the observed 
thermal and scattered light ( §4.11) . 

8. Determine that a tw o-zone model is consistent with steady-state grain dynamics dom- 
inated by collisions kcod type B; 

9. Show that, by relying on models of the inner disk (r < 35 AU) structure, we cannot 
place strong limits on the geometric vertical optical depth of small (0.05 /^m) porous 
grains interior to the birth ring and therefore cannot yet constrain from this ap- 
proach alone ( §4.21) . and 

10. Show that the blow-out size for porous grains is consistent with modeled grain size 
distribution for < lO^M© {UM- 



The next modeling steps will seek to utilize all available imaging and SED data. We 
note in particular that measurements of scattering in polarized light provide strong, comple- 
mentary constraints on the composition and distribution of dust in this system. Future work 
will incorporate fitting to scattered light images rather than surface brightness profiles, and 
will begin to probe the vertical structure of the disk while taking into account the effects 
of projection and blurring by the PSF. Finally, we emphasize that the global structure of 
the debris can be explained by steady-state dust production and diffusion, and we need not 
invoke a planet to clear the inner region of dust. However, the mechanisms responsible for 
maintaining the structure of parent bodies in the birth ring as well as the dust substructure 
remain undetermined. 
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A. ROLL SUBTRACTION FOR AO 

A decisive arbiter of subtraction fidelity is the stability of the stellar PSF. We refined the 
roll subtraction process (Fig. [T]) to simply accommodate the time variability of ground-based 
PSFs. In general, the long-exposure AO PSF is a non-trivial function of atmospheric turbu- 
lence fiuctuations, filtration by closed- loop AO, and additional quasi-static wavefront errors 
arising from the optics in the system ("static speckles"). With the PSF structure arising 
from these effects in mind, we built an empirical model of the PSF for each exposure, and 
then used these models to optimally subtract the stellar image. Our proposed deconstruction 
hinges on the removal of the radial profile from each image. It is much more challenging to 
remove the stellar profile from images with face-on disks. The procedure outlined here is 
most-suited for imaging edge-on disks from the ground. 



A.l. PSF Model 

The ideal model PSF will reproduce the response of the imaging system to the starlight 
which has been occulted by the coronagraph. For our subtraction problem, we are interested 
in estimating the structure of the stellar PSF's outer regions — outside of any coronagraphic 
focal-plane mask. For the moment, we ignore the focal-plane mask and consider only the 
case of a monochromatic image of an on-axis star. We approximate the adaptively corrected 
long-exposure point spread function S'(x) as the convolution {-k) betwee n a "static" PSF a. nd 



a blurring kernel encompassing the time-variable wavefront errors (e.g. IVeran et al.lll997l ) 



S(x)^5s(x)*A(x). (Al) 

Here, A(x) is determined by the variable wavefront errors of partially corrected atmospheric 
turbulence which, depending on correction level, has a characteristic width ^ A/rg (where 
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vq is the Fried length). S's(x) is the PSF that would arise solely from time- independent 
wavefront errors. For a circular pupil of diameter D with no obscurations or wavefront 
errors, 5's(x) is the diffraction-limited Airy pattern of characteristic width ~ X/D. 

With good AO correction, the structure of the long-exposure PSF is characterized by a 
diffraction-limited core and a seeing-disk halo. Features like diffraction spikes and individual 
speckles also play an important role in the structure. To further motivate our model, we 
first decompose the static PSF into the sum of its Airy core {6 < 1.22X/D) and all other 
static speckles {9 > 1.22A/D), 

^s(x) = S's,e(x) + S's,sp(x). (A2) 

Note that since the core is sharply peaked relative to the halo scale in the blurring kernel 
{\/D <^ A/ro), 5's,c(x) A(x) is approximately proportional to A itself; the scale factor a 
is related to the Strehl ratio. For computational simplicity, we treat the blurring effect of 
A(x) on S's,sp(x) by a scalar factor (3, which also correlates with correction level — decreased 
blurring increases peak speckle intensity. Further approximating (]A1|) . we have 

5(x)~aA(x)+/55s,sp(x). (A3) 



We have simplified the PSF from a convolution of static and atmospheric terms (Eq. lAll) 
to a linear combination (Eq. IA3p . and we now seek to increase computational efficiency 
further. First, we note that with this linear combination, it is straightforward to generalize 
to wide-band imaging since A and Sg^sp can each be represented by a linear combination of 
quasi-monochromatic PSFs. Second, although not strictly true under real conditions, we 
further assume that the blurring is symmetric in azimuth (0). This leads us to isolate the 
azimuthally symmetric radial profile of each term in equation (lA3p . For an exposure i, we 
estimate the image of the star as 

^,(^,0)^p,(^)+A5s,sp(^,0), (A4) 

with S's .,p(x) given by S's,sp(x) with the radial profile removed and subsumed into the profile 
term p(x). 

We treat the coronagraphic focal mask in our model (Eq. IA4|) simply. Prima facie, 
the circular mask's suppression of the stellar image will be absorbed into p. Errors in 
centering the st ar behind the mask will result iri additional leakage of light outside the spot 



via diffraction (ILloyd fc Sivaramakrishnarul2005l ). an effect we ignore. These misalignments 



manifest themselves as subtraction errors in the vicinity of the mask edge. 

A key requirement for our model is the profile- removed speckle map, 5*^ sp(x)- practice, 
we use the target images themselves to construct S'sgp(x). We remove the radial profile from 
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each image, mask out the region around the (rotating) edge-on disk, and combine the results 
to form the speckle map. We additionally mask the the diffraction spikes from each image 
when estimating the radial profiles p,, because the spikes' blurring (which we treat only by a 
scale factor /?) can bias the profile estimates. In practice, we mask the disk with a O'.'5-wide 
strip, and the spikes with regions that are 1'.'6 wide, tapering inward for 9 < 5". With Pi(x) 
and gp(x) in hand, the PSF model for each image is parametrized solely by which scales 
the speckle map. 

A. 2. Subtraction Procedure 

We subtract the model PSF Si from each image. The map scale parameter (3i and 
a registration offset arc tuned to minimize subtraction residuals, which are rotated to a 
common frame and combined to form the image of the disk. 

Special software masks are used in various stages of the subtraction process. Diffraction 
spikes and the edge-on disk were masked during calculation of the radial profile. The disk 
was also masked from individual frames before they were combined to form the speckle map. 
Diffraction spikes, the focal plane spot, and the disk were excluded from consideration during 
the least-squares optimization of the PSF model subtraction. Finally, diffraction spikes were 
masked from the subtraction residuals before they were combined to form the final disk 
image. 

In summary, the sequence is to 

1. register each frame to a fiducial position using stellar centroid estimates, 

2. mask the diffraction spikes and the disk, then extract the radial profile about the 
fiducial point (pj), 

3. remove the profile and combine the residuals from each image to obtain the speckle 
map (-S^^gp), 

4. construct a PSF model for each image, using a scaled speckle map {PiS'^^^^) and the 
radial profile (pi), 

5. use an optimization method to subtract PSF model from image, solving for stellar 
position and speckle map scale {(3i), 

6. repeat the process from step (1) using updated centroid estimates if not in final itera- 
tion. 
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7. mask diffraction spikes, rotate the subtraction residuals into the sky's frame, and com- 
bine to form the disk image estimate. 

We find that, for these data, three iterations of this process is enough for convergence to 
visually detectable levels. 



The subtraction procedure described in the previous sections is compared with similar 
implementations in Figure fT2l In panel (a), we follow the steps we outlined above. In panel 
(b), we restrict the model fitting process so that the scale of the speckle map is fixed and 
only offsets were optimized. Finally, in (c) we discard the PSF models altogether and fit the 
offset and scale of an average PSF. No spike or object masking were used when combining 
residual images in (c). In panel (d) we display curves of annular rms (corresponding to the 
photometric errors shown in Figure , excluding the regions near the disk. These curves 
show that fitting the speckle scale gives modest contrast gain (< 0.1 mag), while larger 
gains are obtained when removing the radial profile of the observations prior to combination 
(~ 0.5 mag). 



When quoting sensitivity limits for the detection of point sources around a star, it is 
common practice to measure the standard deviation of pixel values in concentric annuli, cT(r), 
and to set a scale to this curve with artificially inserted sources. In this case, a represents a 
by-eye threshold. Here we wish to extend this methodology for the case of sources within the 
disk. To this end, we adopt a model for the sensitivity and solve for this model's parameters 
through visual detection of randomly inserted sources (both inside and outside of the disk) 
and maximum likelihood techniques. 

We assume that the process of detecting a source can be characterized by zero-mean, 
normally distributed background noise fluctuations. In this case, the background can be 
characterized by a, and the probability of detecting a source (denoting D = 1 for detection) 
of flux / is given by 



which is the cumulative normal distribution function. Our task is to determine a, which 
we can assume is proportional to cr{r) in regions devoid of disk light. Similarly, in the disk 



A. 3. Technique Comparison 



B. MIDPLANE DETECTION SENSITIVITY 
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midplane, we expect an additional contribution to the background fluctuations from the disk 
light. For our ground-based observations, these fluctuations are dominated by speckle noise 
(including Strehl fluctuations) rather than photon statistics, and therefore the variance of this 



contr i bution is proportional to th e square of the disk brightness, P (e.g. lAime fc Soummer 



2004 : iFitzgerald fc Graham! |2006| ). We model the sensitivity in the midplane by 



a^(x,/.,0 = /^'^'(x)+^^/^(x). (B2) 

We fix ^ = to model the sensitivity in regions devoid of disk light. We measure o"(r) with 
a sigma-clipped sample standard deviation of pixels in concentric annuli, excluding regions 
of disk emission. 

We randomly generated positions and fluxes for artificial sources to be inserted in each of 
our ground-based images. Two populations of sources were generated — one set randomly 
distributed along the disk midplane, and another in the off-disk region. The number of 
sources in each population was drawn from a Poisson distribution, and no blending of sources 
was allowed. Fluxes were drawn from a log-uniform distribution about cr(x, /i = 1,^ = 1) in 
the midplane and cr(x, /i = 1,^ = 0) elsewhere. After the computer inserted the randomly 
generated sources, the images were inspected and detections recorded. False detections were 
ignored in our analysis. 

With these data in hand, the problem reduces to finding the most likely model param- 
eters (yU, ^) given the set of artificial source positions, fluxes, and detections {xj,/j,Dj}. 
Our sensitivity curves will use the model parameters which maximize the probability density 
p(yU, ^1 {xj, fi,Di}). This quantity can be re- written with Bayes' Theorem, 

P{f^, ^\ {^iJh Di}) oc p(/i,Op({xi,/i, A} (B3) 

i 

p{fi,OllPiD,\^iJi,^^,0- (B5) 

i 

The final step is valid since we are free to choose the position and flux distribution of 
independent artificial sources without regard to the model parameters. Our prior information 
on the model parameters is represented by p{^,^), and we assume uniform distributions which 
are also independent of /x, ^. We maximized the logarithm of this function with respect to 
/i and ^ for each of our JHK' images, to obtain best-fit l-cr sensitivities, and we show our 
5-0" point source detection sensitivity as a function of separation from the star in Figure O 
Midplane sensit ivities are aver a ged o ver both ansae. We show predicted planet brightnesses 



using models of iBurrows et al.l (119971 ) at different ages, allowing inference of detection limits 
for planet mass. 



-37- 



REFERENCES 

Aime, C. & Soummer, R. 2004, ApJ, 612, L85 

Artymowicz, P. 1988, ApJ, 335, L79 

Augereau, J. C. & Beust, H. 2006, A&A, 455, 987 

Augereau, J. C, Nelson, R. R, Lagrange, A. M., Papaloizou, J. C. B., & Mouillet, D. 2001, 
A&A, 370, 447 

Backman, D. E. & Paresce, F. 1993, in Protostars and Planets III, 1253-1304 

Barrado y Navascues, D., Stauffer, J. R., Song, I., & Caillault, J.-P. 1999, ApJ, 520, L123 

Bryden, G., et al. 2006, ApJ, 636, 1098 

Burrows, A., et al. 1997, ApJ, 491, 856 

Chen, C. H., et al. 2005, ApJ, 634, 1372 

Chen, C. H., et al. 2006, ApJS, 166, 351 

Cutispoto, G., Messina, S., & Rodono, M. 2003, A&A, 400, 659 

Decin, G., Dominik, C, Waters, L. B. F. M., & Waelkens, C. 2003, ApJ, 598, 636 

Dohnanyi, J. S. 1969, J. Geophys. Res., 74, 2431 

Draine, B. T. & Li, A. 2001, ApJ, 551, 807 

Elias, J. H., Frogel, J. A., Matthews, K., & Neugebauer, G. 1982, AJ, 87, 1029 
Fitzgerald, M. P. & Graham, J. R. 2006, ApJ, 637, 541 

Fortney, J. J., Marley, M. S., Hubickyj, O., Bodenheimer, P., & Lissauer, J. J. 2005, As- 
tronomische Nachrichten, 326, 925 

Ghese, W. & Jahreiss, H. 1995, VizieR Onhne Data Catalog, 5070, 

Golimowski, D. A., et al. 2006, AJ, 131, 3109 

Graham, J. R., Kalas, P. G., & Matthews, B. C. 2007, ApJ, 654, 595 
Grigorieva, A., Artymowicz, P., & Thebault, P. 2007, A&A, 461, 537 
Haisch, Jr., K. E., Lada, E. A., & Lada, C. J. 2001, ApJ, 553, L153 



-38- 

Hauschildt, P. H., AUard, F., & Baron, E. 1999, ApJ, 512, 377 
Joy, A. H. & Abt, H. A. 1974, ApJS, 28, 1 
Kalas, P. 2005, ApJ, 635, L169 

Kalas, P., Fitzgerald, M. P., & Graham, J. R. 2007, ApJ, 661, L85 

Kalas, P., Graham, J. R., Clampin, M. C, & Fitzgerald, M. P. 2006, ApJ, 637, L57 

Kalas, P., Liu, M. C, & Matthews, B. C. 2004, Science, 303, 1990 

Keenan, P. C. 1983, Bulletin d'Information du Centre de Donnees Stellaires, 24, 19 

Kenyon, S. J. & Bromley, B. C. 2004, AJ, 127, 513 

Krisciunas, K., et al. 1987, PASP, 99, 887 

Krist, J. E., et al. 2005, AJ, 129, 1008 

Krist, J. E. & Hook, R. 2004, The Tiny Tim User's Guide, Version 6.3 (Baltimore: STScI) 
Kuchner, M. J. & Holman, M. J. 2003, ApJ, 588, 1110 

Lagrange, A.-M., Backman, D. E., & Artymowicz, P. 2000, Protostars and Planets IV, 639 
Lecaveher Des Etangs, A., Vidal-Madjar, A., & Ferlet, R. 1996, A&A, 307, 542 
Li, A. & Greenberg, J. M. 1998, A&A, 331, 291 

Linsky, J. L., Bornmann, P. L., Carpenter, K. G., Hege, E. K., Wing, R. F., Giampapa, 
M. S., & Worden, S. P. 1982, ApJ, 260, 670 

Liou, J.-C. & Zook, H. A. 1999, AJ, 118, 580 

Liu, M. C. 2004, Science, 305, 1442 

Liu, M. C, Matthews, B. C, Williams, J. P., & Kalas, P. G. 2004, ApJ, 608, 526 
Lloyd, J. P. & Sivaramakrishnan, A. 2005, ApJ, 621, 1153 

Marley, M. S., Fortney, J. J., Hubickyj, O., Bodenheimer, P., & Lissauer, J. J. 2007, ApJ, 
655, 541 

Marois, C, Lafreniere, D., Doyon, R., Macintosh, B., & Nadeau, D. 2006, ApJ, 641, 556 



-39- 

Masciadri, E., Mundt, R., Henning, T., Alvarez, C, & Barrado y Navascues, D. 2005, ApJ, 
625, 1004 

Mathioudakis, M. & Doyle, J. G. 1991, A&A, 244, 433 
Mathis, J. S. & Whiffen, G. 1989, ApJ, 341, 808 

Metchev, S. A., Eisner, J. A., Hillenbrand, L. A., & Wolf, S. 2005, ApJ, 622, 451 

Meyer, M. R., Backman, D. E., Weinberger, A. J., & Wyatt, M. C. 2007, in Protostars and 
Planets V, ed. B. Reipurth, D. Jewitt, & K. Keil, 573-588 

Moro-Martm, A. & Malhotra, R. 2005, ApJ, 633, 1150 

Moro-Martin, A., Wolf, S., & Malhotra, R. 2005, ApJ, 621, 1079 

Neuhauser, R., Guenther, E. W., Alves, J., Huelamo, N., Ott, T., & Eckart, A. 2003, 
Astronomische Nachrichten, 324, 535 

Okamoto, Y. K., et al. 2004, Nature, 431, 660 

Ozernoy, L. M., Gorkavyi, N. N., Mather, J. C, & Taidakova, T. A. 2000, ApJ, 537, L147 
Pantin, E., Lagage, P. O., & Artymowicz, P. 1997, A&A, 327, 1123 
Perryman, M. A. C, et al. 1997, A&A, 323, L49 

Persson, S. E., Murphy, D. C, Krzeminski, W., Roth, M., & Rieke, M. J. 1998, AJ, 116, 
2475 

Pinte, C, Menard, P., Duchene, G., & Bastien, P. 2006, A&A, 459, 797 
Plavchan, P., Jura, M., & Lipscy, S. J. 2005, ApJ, 631, 1161 

Roberge, A., Weinberger, A. J., Redfield, S., & Feldman, P. D. 2005, ApJ, 626, L105 
Robinson, R. D., Linsky, J. L., Woodgate, B. E., & Timothy, J. G. 2001, ApJ, 554, 368 
Roques, F., SchoU, H., Sicardy, B., & Smith, B. A. 1994, Icarus, 108, 37 
Sheehy, C. D., McCrady, N., & Graham, J. R. 2006, ApJ, 647, 1517 
Skrutskie, M. P., et al. 2006, AJ, 131, 1163 

Song, I., Weinberger, A. J., Becklin, E. E., Zuckerman, B., & Chen, C. 2002, AJ, 124, 514 



-40- 

Strubbe, L. E. & Chiang, E. I. 2006, ApJ, 648, 652 

Telesco, C. M., et al. 2005, Nature, 433, 133 

Thommes, E. W. & Lissauer, J. J. 2003, ApJ, 597, 566 

Torres, C. A. O. & Ferraz Mello, S. 1973, A&A, 27, 231 

Veran, J.-R, Rigaut, F., Maitre, H., & Rouan, D. 1997, 14, 3057 

Voshchinnikov, N. V., H'in, V. B., Henning, T., & Dubkova, D. N. 2006, A&A, 445, 167 
Voshchinnikov, N. V. & Mathis, J. S. 1999, ApJ, 526, 257 

Wahhaj, Z., Koerner, D. W., Ressler, M. E., Werner, M. W., Backman, D. E., & Sargent, 
A. I. 2003, ApJ, 584, L27 

Weinberger, A. J., Bccklin, E. E., & Zuckcrman, B. 2003, ApJ, 584, L33 

Wyatt, M. C. 2006, ApJ, 639, 1153 

Wyatt, M. C, Dermott, S. F., Telesco, C. M., Fisher, R. S., Grogan, K., Holmes, E. K., & 
Piiia, R. K. 1999, ApJ, 527, 918 

Zuckerman, B. 2001, ARA&A, 39, 549 

Zuckerman, B., Forveille, T., & Kastner, J. H. 1995, Nature, 373, 494 
Zuckerman, B., Song, I., Bessell, M. S., & Webb, R. A. 2001, ApJ, 562, L87 



This preprint was prepared with the A AS IM^jX macros v5.2. 



-41 - 



Table 1. Observations of AU Mic 



night 


band 


'"mask ( ) 


r(s) 


APA (°) 


2004 Aug 29 


J 


0.75 


600 


16.9 




H 


0.75 


690 


49.8 




K' 


0.75 


1110 


40.5 


2004 Aug 30 


J 


0.50 


570 


39.0 




H 


0.50 


600 


34.5 




K' 


0.50 


1260 


45.6 




n 


0.50 


720 


43.6 



Note. — The radius of the coronagraphic focal 
plane mask is given by rmask- Filters were cycled after 
a few short exposures in each band. Here, T is the 
total integration time and APA is the total amount of 
field rotation over the course of the exposure sequence 
in that band. 
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Table 2. Midplane Surface Brightness Power-Law Indices 



band ansa fit domain 

15-32 AU 32-60^ AU 



F606W 


NW 


1, 


.46 


± 


0. 


,09 


3, 


.00 


± 


0, 


,13 




SE 


1, 


.53 


± 


0. 


,08 


2, 


.97 


± 


0, 


,12 




avg. 


1, 


.49 


± 


0. 


,09 


2, 


.99 


± 


0, 


,12 


J 


NW 


1, 


.21 


± 


0. 


,10 


3, 


.95 


± 


0. 


,19 




SE 


1, 


.34 


± 


0. 


,12 


3, 


.27 


± 


0. 


,19 




avg. 


1, 


.27 


± 


0. 


,11 


3, 


.61 


± 


0, 


,19 


H 


NW 


1, 


.19 


± 


0. 


,19 


3, 


.84 


± 


0, 


,13 




SE 


1, 


.58 


± 


0. 


,15 


3, 


.47 


± 


0, 


,17 




avg. 


1, 


.39 


± 


0. 


,17 


3, 


.66 


± 


0. 


,15 


K' 


NW 


1, 


.24 


± 


0, 


,13 


5, 


.16 


± 


0, 


,23 




SE 


1, 


.09 


± 


0. 


,10 


4, 


.63 


± 


0, 


,16 




avg. 


1, 


.17 


± 


0. 


,12 


4, 


.90 


± 


0, 


,20 



1'Maximum projected distance used in fit; 
the domain may be further restricted by the 
availabihty of data in Fig. [51 

Note. — Power-law indices a, calcu- 
lated by converting the the midplane surface 
brightness profiles (Fig. [5]) to flux units and 
fitting f[h) oc 6^". Formal 1-a errors are 
scaled by when this quantity is > 1. 

These errors are lower limits, as both sys- 
tematic errors and measurement covariance 
have been ignored. Entries marked "avg." 
are computed by averaging a and its vari- 
ance over both ansae. 
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Table 3. Comparison of Disk Features 



Label 


NW 


SE 




Location (AU) 








L04 


K05 


M05 this work 


A 


T 


T 


25 


26 


22 25 


B 




i 


28.5 


29 


26 28 


C 




T 


31 


33 


32 32 


D 


t 


i 


37 


37 


38 37 


E 


i 


i 




23 


21 



Note. — Features are marked in Figure [71 Ar- 
rows denote localized enhancements (|) an d defi cits 



(I) of disk brightness. Featu re lo catio ns in iKOSi are 
seen in F606W, while |L04| and iMOSi are H-hsnad. 
Average feature locations in the near-IR data from 
this paper are in the final column, based on visual 
positions in Fig. [7] with uncertainties of approxi- 
mately ±0.5 AU. We newly designate the close-in 
feature E here. 
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Table 4. Best-Fit Model Parameters 



region range (AU) dust mass (M^) a IP 



porous fMathis fc Whiffen 19891 



inner 
outer 



35-40 
40-300 



compact fPraine fc Li 2001') 
inner 28-32 
outer 32-300 



1.0x10-2 
2.3x10"^ 



2.7x10^2 
2.7x10"^ 



3-6 mm -3.5^ +1.5 

0.05-3.0 /im -4.1 -2.5 

1-2 mm -3.5^ +1.5 

0.15-50.0 /im -4.1 -2.5 



*Fixed. 



Note. — Parameters for our models which produce the best fit to the F606W- 
and i/-band profiles and SED, for both porous and compact grains. The grain size 
distribution is dN{a) oc a''da, while the surface density is described by S(r) oc r'^. 
Illuminating star: T^q =3600 K, log (7 = 4.5 NextGen model, stellar radius 
0.88 i?0. Both models give comparable fits to the scattered light and SED, but 
the measured polarization rules out compact /im-sized grains (Figs. fTOl fc fTTi) . 
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Table 5. Best-Fit Model Avg. Scattering Parameters 



A g A 



porous 

606 nm 0.83 0.52 

1.6/im 0.81 0.54 

compact 

606 nm 0.66 0.83 

1.6/im 0.60 0.81 



Note. — Scatter- 
ing parameters for the 
outer regions of our 
best-fit models (Ta- 
ble H]), averaged over 
the grain size distri- 
bution. The scatter- 
ing asymmetry factor 
g and albedo A are 
given. 
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Fig. 1. — The roll subtraction process, which estimates the image of the disk in the presence 
of a fluctuating stellar PSF. The disk (a) is observed on an alt-az telescope with no image 
de-rotation (b). The PSF is fixed relative to the detector while the disk image rotates due 
to sidereal motion. The observed disk images are blurred by different amounts by the time- 
varying PSF, illustrated by different line weights the disk images and the size of the circular 
stellar halo. Estimates of the PSF scaled and positioned to match the stellar image (c) 
are subtracted from the observations, leaving residual disk images (d). These images are 
transformed to a common frame (e) and combined to estimate the disk image (f). 



60 



40 



20 







20 



40 



60 



AU 




10- 



10-4 



10-^ 



I /F^ (arcsec ^) 



Fig. 2. — Images of the AU Mic debris disk with direct starhght removed via PSF roll 
subtraction ( §2.2p . Data are displayed on a square- root scale, and in each band the disk 
brightness I is divided by the stellar flux F^, allowing for comparison of color differences 
intrinsic to the dust rather than the star. In each image, a circular software mask (I'.'S 
diameter) is apphed to obscure subtraction residuals about the stellar location, marked by 
a small white circle. The near-IR images have been additionally smoothed by a Gaussian 
matched to the resolution of the PSF to eliminate small high-frequency errors introduced 
by the masking process (Appendix El) . The data in panel (a) were obtained with the ACS 
coronagraph aboard HST, while the data in (h)-(d) are newly acquired via Keck AO. The 
blue color of the dust is indicated by the trend of decreasing scattering efficiency toward 
longer wavelengths in panels (a)-(d). 



-48- 




n 



Fig. 3. — A JHK' composite image using the data in Figure O The data are displayed 
relative to the stellar brightness, highlighting the intrinsic color of the dust. The structures 
emanating from the mask outside of the disk are residuals from the stellar PSF subtraction. 
The blue color of the disk is visible, as is evidence for substructure. 
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Fig. 4. — Vertical structure in the near-IR disk images. Panels (a) and (b) show a JHK' 
composite of disk images, whereby a Gaussian fit to the vertical structure of the disk has 
been subtracted. The residuals clearly show the location of the sharp midplane, and have 
been scaled by a smooth fit to the amplitudes of the Gaussian functions. The SE side is 
flipped about the star to allow for direct comparison to the NW. The locations of physical 
features are indicated (A-E; cf. Table [3]). Panels (c) and (d) plot the variation of the 
vertical midplane position resulting from Gaussian fits to the vertical profile. No significant 
differences in midplane position with wavelength are seen. 
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Fig. 5. — Midplane surface brightness profiles for the disk ansae. Flux was gathered in 
O'.'l X O'.'S apertures, which were placed according to the spline fit to the vertical disk mid- 
plane. The measurements are shown relative to the stellar brightness to highlight the intrinsic 
scattering properties of the dust in the J- to K'-hands. The uncertainties represent 1-a ran- 
dom errors, and do not include systematic errors in calibration ( §2.31) . We also indicate 3-a 
upper limits to the L' brightness. 
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Fig. 6. — The disk color vs. the projected distance along the disk midplane. We compute 
disk colors relative to the i/-band midplane surface brightness profile. The contribution of 
the intrinsic stellar color has been removed. The disk appears blue between the visible and 
near-IR and between H and K' . Note that K'-H colors are shown, which accounts for the 
reversed gradient from F606W-H. Systematic uncertainties in calibration are not included, 
though these are expected to be < 0.3 mag for K'-H and J-H, and < 0.2 mag for F606W- 
H ( §2.3f[27il) . In the inner disk (< 35 AU), the colors are consistent with a fiat profile. The 
outward branching of these curves in the outer disk indicates blue color gradients (most 
clearly indicated by the F606W-K' color). These gradients may indicate differences in grain 
sizes or compositions. 
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Fig. 7. — Substructure in the disk images. The upper and lower panels are F606W data, 
while the middle panels are JHK' composites. The image in each band is scaled by a spline 
fit to an average of its NW and SE midplane surface brightness profiles ( §3.2p . The SE 
side is flipped about the star to allow for direct comparison to the NW. Feature locations 
in the near-IR are indicated by so lid lines (A-E), while the original identification of feature 
locations in F606W data by lK05l are shown as do tted lines (cf. Table [3]) • A broad clump in 
the SE ansa at 48 AU was also identified by lK05l . We confirm the presence of a brightness 
deficit at location E. In our favored models (|SD, the bulk of the scattered light comes from 
small grains outside of 40 AU. Therefore, the features at projected distances inward of 40 AU 
(A-E) must arise from azimuthal perturbations in the dust distribution outside this radius. 
The origin of the substructures is unknown; they may result from the gravitational influence 
of unseen planets. 
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Fig. 8. — Substructure in disk profiles. Surface briglitness profiles computed with C'l x O'.'l 
apertures have been processed in the same manner as Fig. [TJ by dividing each wavelength's 
profile by a smooth spline function derived from an average of the NW and SE profiles from 
Fig. [51 Feature locations derived for the near-IR data are indicated by s olid v ertical lines 
(A-E), while dotted lines indicate the feature locations in F606W data by IK05 . 
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Fig. 9. — Point source detection limits (5-cr). The solid curves are for sources outside 
of the disk, while the dashed curves give limits for sources residing in the disk midplane. 
In some places the midplane is more than a magnitude less sensitive. No point sources 
were detected along t he disk. Horiz o ntal l ines indicate the predicted brightnesses for the 
model giant planets of Isurrows et aL ( 1997 ). at ages of 10 and 30 Myr (AU Mic is 12^4 Myr; 
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Fig. 10. — A model fit to the scattered light and SED of the system, (left) The surface 
brightness profiles from Fig. O along with surface brightness profiles from the best-fit model. 
The gray boxes above the profiles represent the grain locations in our model (cf. Table Hj); 
the dark region indicates the inner region of larger grains, while the smaller scatterers are 
in the lighter zone outside, (right) The model SED along with measured photometry of 
AU Mic. In this model, the smaller grains are responsible for the bulk of the scattered light 
and the mid-IR emission, while the larger grains reproduce the long-wavelength end of the 
SED. 
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Fig. 11. — (left) A comparison of the surface brightness profiles from Fig. [5] to surface 
brightness profiles from the best-fit model in J and K' bands. These data were not used 
in the fitting process. The differences in overall flux may be due to photometric calibration 
uncertainties ( §2.31) . though the model clearly overestimates the emission in the outer zone in 
these bands. More complicated models may resolve these discrepancies; a changing minimum 
grain size with radius can produce a color gradient, (right) The fra ction of linear F606 W 
polarization produced by our model compared to measurements of iGraham et al.l (120071 ). 
which include both systematic and random errors. Compact, /im-sized grains are ruled out 
by these data. 




.tl I - 





CaJ 




m 


- V \ 


- -Cc) 






1.1. 





1 2 :5 4 5 
aeparation C'') 



Fig. 12. — A comparison of different methods of roll subtraction. Panel (a) shows the 
resulting image from the technique described in §2.2[ Panel (b) is the same as (a), except 
the scale of the speckle map is not optimized ( §A.2p . For comparison, panel (c) registers 
an average PSF to each image, rather than fitting a PSF model. Each of these images is 8" 
on a side. In (d), we show the annular rms (excluding the disk) of photometry in O'.'l x 0'.'5 
apertures as a function of radius for subtractions shown in (a)-(c). At a radius of 2", curve 
(b) is 0.56 mag more sensitive than (c), while (a) is 0.16 mag more sensitive than (b). The 
gain in contrast when subtracting the profile highlights the suitability of the technique for 
edge-on disks. 



